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This  dissertation  describes  the  application  of  endo- 
chronic  ( endos=internal , chronos=time ) theory,  which  is  a 
unified  viscoplastic  theory,  to  the  analysis  of  fiber 
reinforced  concrete  structures  subjected  to  marine  exposure 
Endochronic  analysis  has  opened  up  new  avenues  of 
constitutive  modeling  enabling  a realistic  analysis  of  struc 
tures  taking  into  account  material  nonlinearities.  Unlike 
the  conventional  plasticity-based  theories,  this  formulation 
does  not  need  the  definition  of  a yield  surface;  instead  it 
uses  the  concept  of  an  intrinsic  time  defined  in  terms  of 
both  strain  and  time  to  analyze  elastic  as  well  as  plastic 
behavior  in  a continuous  manner.  The  theory  has  been 
successfully  used  to  analyze  metals,  soils,  and  concrete. 

The  experimental  data  are  obtained  by  testing  fiber 
reinforced  concrete  specimens  subjected  to  varying  periods 
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of  circulating  sea  water  exposure.  Uniaxial  compressive 
stress-strain  behavior  is  determined  by  testing  concrete 
cylindrical  specimens  in  parallel  with  a co-axial  special 
alloy  heat  treated  steel  tube.  The  endochronic  constants 
defining  the  material  behavior  are  found  by  using  an  opti- 
mization procedure  based  on  the  experimental  stress-strain 
behavior.  Experimentally  determined  flexural  behavior  is 
predicted  based  on  the  endochronic  constants  and  a proba- 
bilistic method  to  account  for  the  tensile  contribution  of 
fibers  intersecting  a unit  cross  sectional  area.  The  agree- 
ment between  the  experimental  and  predicted  behavior  is 
very  encouraging. 

The  presentation  includes  a description  of  the  beha- 
vior of  concrete  in  the  marine  environment,  the  properties 
of  fiber  reinforced  concrete,  and  recent  improvements  in  the 
endochronic  formulation,  such  as  plastic  strain  and  jump 
kinematic  hardening  formulations,  which  overcome  the  earlier 
deficiencies  like  the  open  hysteresis  loops  under  cyclic 
loading.  A simplified  extension  of  the  endochronic  formu- 
lation is  presented  to  analyze  the  thermally  sensitive  mate- 
rials like  ice  and  frozen  soil.  This  formulation  is  used 
in  conjunction  with  a finite  element  analysis  program  to 
analyze  the  response  of  a typical  Arctic  offshore  island 
foundation.  The  numerical  procedures  used  are  described  in 
the  Appendix. 
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CHAPTER  I 
INTRODUCTION 


1 . 1 General  Introduction 

This  dissertation  is  an  analytical  and  experimental 
study  of  fiber  reinforced  concrete  exposed  to  the  marine 
environment  and  frozen  soil.  An  endochronic  (endos=  inter- 
nal, chronos=time ) theory,  a unified  viscoplastic  theory, 
is  used  to  analytically  predict  the  experimental  flexural 
behavior.  The  endochronic  formulation  has  been  extended  by 
including  the  thermal  and  creep  effects  to  analyze  the 
response  of  a typical  Arctic  offshore  island  foundation. 

Concrete  is  generally  considered  to  be  an  ideal 
marine  construction  material,  requiring  little  or  no  main- 
tenance, due  to  its  inertness  in  sea  water.  Good  concrete 
with  a low  permeability,  submerged  in  sea  water,  increases  in 
strength  with  the  passage  of  time.  Concrete  has  been  used  in 
the  past  for  marine  structures  like  breakwaters  and  sea  walls, 
and  even  in  ship  building.  With  the  increase  in  offshore 
oil  exploration  activity,  concrete  is  assuming  an  increased 
role  in  offshore  construction.  Concrete  gravity  platforms 
are  very  well  adapted  for  the  severe  environmental  areas  like 
the  North  Sea  and  Arctic  regions.  Because  of  the  long  lasting 
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nature  of  concrete  in  sea  water,  nuclear  industry  is  consi- 
dering the  use  of  thick  concrete  shells  for  marine  disposal 
of  nuclear  wastes. 

In  spite  of  its  high  compressive  strength,  concrete 
has  low  tensile  strength.  This  is  the  major  factor  in  the 
initiation  of  cracks  and  deterioration  of  many  reinforced 
concrete  structures  in  the  marine  environment.  When  the 
reinforcing  steel  corrodes,  the  resulting  oxide  products 
occupy  larger  volumes  giving  rise  to  internal  tensile  stre- 
sses. Once  the  concrete  cracks,  the  corrosion  mechanism 
accelerates  due  to  the  easier  access  of  the  corroding  medium 
to  the  embedded  steel.  This  phenomenon  has  lead  to  the 
deterioration  of  many  coastal  structures.  The  cost  of 
replacement  and  repair  of  bridge  structures  alone  in  the 
U.S.A.  is  projected  to  be  around  6.3  billion  dollars.  Some 
of  the  methods  used  to  prevent  the  cracking  from  the  embe- 
dded metal  corrosion  are  i)  use  of  high  strength  concrete 
with  larger  clear  cover,  ii)  use  of  epoxy  coated  steel  reinfor 
cement,  iii)  chemical  inhibitors,  and  iv)  cathodic  protection. 
High  strength  concrete  has  a lesser  porosity  affording  a great 
er  immunity  to  the  steel  from  the  corroding  medium.  Research 
is  being  done  on  the  use  of  epoxy  coatings  for  steel,  and  some 
corrosion  inhibitors  like  calcium  nitrate  in  the  concrete 
mix  for  retarding  or  preventing  steel  corrosion. 
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1 .2  Special  Concretes 

Recently,  many  techniques  have  been  developed  to 
improve  the  tensile  behavior  of  concrete.  One  effective  but 
expensive  method  is  to  reduce  the  internal  flaws  by  filling 
the  pores  with  polymeric  substances.  There  are  three  differ- 
ent types  of  polymer  concretes: 

1)  Polymer-Impregnated  Concrete  (PIC):  The  capillary  pores 

in  the  hardened  concrete  are  impregnated  by  a liquid  monomer 
which  is  subsequently  polymerized  to  form  the  solid  polymer 
insitu  within  the  pores. 

2)  Polymer  Concrete  (PC):  Consists  of  a mix  of  aggregate  and 

polymer  with  no  cement.  The  polymer  acts  as  the  binder  for 
the  aggregate. 

3)  Latex-Modified  Concrete  ( LMC ) : Polymer  latex,  a colloidal 

suspension  of  polymer  in  water,  is  added  to  the  concrete  mix 
at  the  mixing  stage. 

All  the  methods  produce  hard,  dense,  and  high 
strength  concrete.  Some  of  the  applications  of  polymer 
concretes  have  been  in  the  overlays  for  pavements  and  in  the 
construction  of  seafloor  spherical  hull  structures  for  oil 
storage . 

Another  way  of  improving  the  tensile  behavior  of 
concrete  is  by  adding  short,  discontinuous  fibers  which 
restrict  the  crack  growth  by  bridging  across  them.  The 
fibers  make  the  concrete  behave  like  a pseudo-ductile  mate- 
rial instead  of  its  usual  brittle  behavior.  They  have  been 
used  to  reinforce  brittle  materials  since  ancient  times; 
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straw  was  used  to  reinforce  sunbaked  bricks  even  during  the 
Biblical  periods,  and  more  recently,  asbestos  fibers  have 
been  used  to  reinforce  Portland  cement.  Fiber  reinforced 
concrete  has  been  applied  in  many  instances  where  the 
performance  of  ordinary  concrete  was  unsatisfactory.  Steel 
fiber  reinforced  concrete  was  found  to  reduce  the  thickness 
of  airport  runways  by  one  half.  It  has  been  used  in  the 
successful  repair  of  stilling  basins  of  dams  damaged  by  float- 
ing debris  digging  through  the  ordinary  concrete  linings. 
Another  successful  application  was  for  the  linings  in  hydraulic 
structures  where  cavitation  was  a big  problem.  Steel  fiber 
reinforced  dolosse  units  were  found  to  survive  severe  winter 
storms  during  which  time  corresponding  ordinary  concrete  units 
suffered  extensive  damage.  Polypropylene  fibers  are  widely 
used  for  concrete  piles,  in  view  of  their  excellent  impact 
resistance  properties.  Glass  fibers  have  been  indicated  as  a 
substitute  for  the  hazardous  asbestos  fibers.  Typical  marine 
application  of  fiber  reinforced  concrete  is  that  for  break- 
water armor  units. 

1 . 3 Constitutive  Modeling  Techniques 

In  the  analysis  of  concrete  structures,  concrete  has 
been  modeled  as  a linear  elastic,  ideally  plastic,  or  work- 
hardening plastic  material.  In  view  of  the  highly  nonlinear 
behavior  of  concrete,  all  the  above  methods  are  only 
approximations.  Linear  elastic  modeling  is  reasonable  only 
at  low  stress  levels;  for  high  stress  levels  plastic  formu- 
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lations  become  necessary.  These  theories,  originally  devel- 
oped for  metals  and  later  adapted  for  concrete,  assume  the 
existence  of  a yield  criterion  coupled  with  a hardening  rule 
to  define  the  subsequent  yield  surfaces.  There  is  no  unique- 
ness associated  with  the  choice  of  yield  surfaces  leading 
to  an  ambiguity  in  the  results.  These  models  can  be  viewed 
as  discontinuous-material  models,  as  the  material  response 
is  separately  analyzed  in  various  stages  i.e.,  loading, 
unloading,  and  reloading.  However,  real  material  behavior 
is  usually  continuous  and  includes  many  complicated  cross 
effects.  Very  often,  such  discontinuous  models,  instead  of 
simplifying  the  problem,  are  the  sources  of  numerical  diffi- 
culties and  inef f iciences . 

On  the  other  hand,  the  recently  developed  endochronic 
theory  describes  the  material  behavior  in  a continuous 
manner  without  the  need  of  a yield  surface  definition.  This 
overcomes  the  uncertainties  associated  with  the  choice  of 
yield  surface.  This  theory  describes  the  behavior  in  terms 
of  an  intrinsic  time  function,  which  is  a non-decreasing 
scalar  variable  that  depends  on  the  increments  of  both 
strain  and  time.  The  intrinsic  time  parameter  represents 
the  extent  of  change  or  damage  of  the  internal  structure  of 
a material  subjected  to  deformation.  This  seems  to  be  one  of 
the  most  comprehensive  theories,  capable  of  modeling  strain 
hardening,  softening,  unloading  and  reloading,  cross  hardening 
and  other  nonlinear  features  of  concrete. 
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1 . 4 Outline  of  the  Dissertation 

Chapter  II  reviews  the  literature  in  the  areas  of 
marine  concrete,  fiber  reinforced  concrete,  and  the  consti- 
tutive theories  to  model  its  behavior.  The  discussion 
includes  the  causes  for  the  deterioration  of  concrete 
structures  in  the  marine  environment,  improvements  in  the 
concrete  properties  due  to  the  fiber  addition,  various  past 
and  present  constitutive  modeling  theories,  and  the  tech- 
niques for  the  analysis  of  concrete  structures.  Chapter  III 
discusses  the  adaptation  of  endochronic  theory  to  model 
various  nonlinear  features  of  concrete.  The  experimental 
techniques  to  obtain  the  complete  stress-strain  curves  and 
the  determination  of  endochronic  constants  are  also  presen- 
ted. The  development  of  an  analytical  model  to  predict  the 
load-deflection  behavior  of  fiber  reinforced  concrete  beams 
is  discussed  in  Chapter  IV.  The  analytically  predicted 
behavior  is  compared  with  that  obtained  by  experimental 
methods.  Chapter  V presents  an  extension  of  the  endochronic 
formulation  to  include  the  thermal  sensitivity  of  materials 
like  ice  and  frozen  soils.  This  formulation  is  used  in 
conjunction  with  a finite  element  analysis  program  to 
analyze  the  foundation  response  of  a typical  Arctic  offshore 
island  composed  of  frozen  soil.  Chapter  VI  concludes  the 
dissertation  with  some  suggestions  for  further  work  in  this 
area.  Various  numerical  techniques  used  in  the  investi- 
gation are  discussed  in  Appendix-A. 


CHAPTER  II 
LITERATURE  REVIEW 


This  chapter  reviews  the  literature  on  four  different 
aspects  covering  the  entire  dissertation  as  follows: 

1 ) Durability  of  concrete  in  sea  water, 

2)  Fiber  reinforced  concrete, 

3)  Constitutive  models  for  concrete, 

and 

4)  Flexural  analysis  of  concrete  beams. 

2 . 1 Durability  of  Concrete  in  Sea  Water 

There  are  basically  three  corrosion  zones  in  a marine 
structure:  immersed,  splash  and  tidal,  and  atmospheric  as 

shown  in  Fig.  2.1.  The  basic  mechanism  of  deterioration  of 
concrete  in  any  zone  is  due  to  an  internal  expansion  of 
either  the  concrete  or  reinforcing  steel.  In  addition,  the 
structural  elements  are  subjected  to  different  environmental 
and  biological  effects  depending  on  the  vertical  location. 

Marine  concrete  structures  can  fail  due  to  the  dete- 
rioration of  the  reinforcing  steel  or  the  concrete  itself. 
Magnesium  salts  dissolved  in  the  sea  water  can  deteriorate 
the  concrete  without  the  corrosion  of  steel  reinforcement. 
Sulphates,  present  in  the  sea  water  mostly  as  magnesium 
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Fig.  2.1  Corrosion  Zones  in  an 

Offshore  Concrete  Structure 


9 


salts,  react  with  tricalcium  aluminate  (C3A)  in  the  cement 
paste  to  form  sulphoaluminate  hydrate  ( ettringite ) , which  has 
more  than  twice  the  volume  of  C3A,  leading  to  a swelling  of 
concrete  and  cracking.  A proper  choice  of  cement  with  low 
C3A  contents,  low  water-cement  ratio,  and  air  entrainment 
reduces  the  damage  due  to  this  kind  of  corrosion.  In  the 
normal  marine  environment  this  effect  is  negligible  com- 
pared to  the  effect  due  to  reinforcement  corrosion. 

When  the  steel  corrodes,  its  corrosion  products  occu- 
py a considerably  larger  volume,  typically  two  to  three 
times  the  original  volume  of  metal  before  corrosion.  As  a 
result,  high  internal  tensile  stresses  are  created  which  can 
disrupt  the  surrounding  concrete  and  cause  cracks  to  form. 
Reduced  cross  sectional  areas  of  steel  may  also  lead  to 
structural  failures,  being  more  critical  with  the  prestress- 
ing steel  tendons  than  with  larger  reinforcing  bars.  The 
alkalinity  of  cement  in  concrete  is  usually  adequate  to  pro- 
vide the  required  pH  level  to  make  the  steel  passive  to 
corrosion.  However,  over  a period  of  time  chlorides  in  the 
sea  water  penetrate  the  concrete  and  depassivate  the  steel, 
making  it  corrode  even  when  the  surrounding  concrete  is 
alkaline.  Another  mechanism  of  steel  corrosion  is  due  to 
the  non-uniformities  in  the  chemical  composition  and  the 
resultant  electric  potentials.  Supply  of  oxygen  is  essen- 
tial for  this  mechanism  to  function.  Houston  and 
Ferguson  (1972)  reported  that  the  ratio  of  clear  cover  to 


bar  diameter  should  exceed  three,  to  provide  a reasonable 
protection  to  steel. 

In  the  immersed  zone  the  structural  members  are  acted 
upon  by  marine  organisms  and  fouling  elements  like  barna- 
cles. They  increase  the  surface  area  and  the  roughness, 
attracting  greater  forces  on  the  structure.  They  could  also 
cause  local  corrosion  and  differential  aeration  cells  due  to 
the  pockets  of  entrapped  oxygen.  A hard,  dense  concrete 
with  a smooth  surface  is  said  to  reduce  the  chances  of  this 
kind  of  a damage. 

In  the  splash  and  tidal  zones,  the  additional  corro- 
sion mechanisms  are  due  to  alternate  wetting  and  drying  from 
tide  cycles,  wear  and  abrasion,  and  a localized  pitting  due 
to  cavitation.  Alternate  wetting  and  drying  may  cause 
cracks  due  to  the  resulting  volume  changes.  Evaporated  sea 
water  usually  leaves  behind  its  salt  contents  which  can 
accelerate  the  corrosion.  Again,  a hard,  dense  concrete  is 
reported  to  reduce  the  damage  in  this  zone. 

Corrosion  in  the  atmospheric  zone  is  mainly  due  to 
the  salt-laden  humid  atmosphere  above  the  sea  surface.  If 
the  concrete  is  permeable  the  chlorides  in  the  atmosphere 
reach  the  steel  reinforcement  in  the  concrete  and  promote 
corrosion.  Another  mechanism,  though  rare  in  marine  envi- 
ronment, can  be  due  to  the  reaction  of  carbon  dioxide  in  the 
atmosphere  with  the  hydroxides  in  the  concrete.  This  reac- 
tion, referred  to  as  carbonation,  reduces  the  alkalinity  of 
concrete  to  such  a level  where  the  passivity  is  destroyed 


and  corrosion  can  start.  Additional  details  on  the  corro- 
sion mechanisms  of  concrete  in  sea  water  have  been  discu- 
ssed by  Beeby  (1978),  Mather  (1982),  Schupak  (1982),  and 
Taylor  and  Sharp  (1978). 

Field  studies  to  determine  the  durability  of  concrete 
have  been  carried  out  by  O'Neil  (1980).  Reinforced  concrete 
beams,  with  and  without  air  entrainment,  were  exposed  to 
natural  weathering  conditions  at  Treat  Island,  Maine,  for 
over  twenty  years.  The  effect  of  various  concrete  mix  para- 
meters on  the  durability  of  concrete  was  evaluated.  Similar 
studies  have  also  been  carried  out  by  Gjorv  (1971)  in  Norway. 
The  nature  of  the  two  investigations  requires  a long  dura- 
tion of  time  to  evaluate  the  durability  of  concrete. 

Nishibayashi , Yamura,  and  Inoue  (1980)  described  an 
accelerated  test  procedure  to  evaluate  the  durability  of 
concrete  in  sea  water.  Concrete  mixes  with  different  types 
of  cement,  water-cement  ratios,  and  curing  histories  were 
subjected  to  cycles  of  immersion  in  sea  water  for  24  hours 
followed  by  an  oven  drying  at  60  to  80°C.  In  addition  to 
visual  observation  for  cracks,  the  reduction  in  dynamic 
modulus  of  elasticity  and  length  changes  were  made  every 
ten  cycles  for  200  cycles.  By  comparing  the  extent  of 
damage  and  other  physical  characterestics  to  the  findings  of 
Gjorv  (1971),  they  concluded  that  200  cycles  of  accelerated 
testing  corresponds  to  a 30-year  exposure  in  the  marine 
environment,  i.e.,  each  cycle  is  equivalent  to  50  to  60 
days  of  actual  field  exposure.  Similar  techniques  were 
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employed  by  Reddy,  Gopal,  Nomikos,  and  Thiel  (1984)  to  eva- 
luate the  durability  of  marine  fiber  reinforced  concrete. 

A numerical  method  to  predict  the  corrosion  rates  of 
reinforcement  and  the  time  required  for  the  first  crack 
formation  was  proposed  by  Bazant  (1979a).  Various  features 
included  in  the  model  are  transport  of  oxygen  and  chloride 
ion  through  the  concrete  cover,  mass  sinks  and  sources 
of  oxygen,  ferrous  hydroxide  and  hydrated  rust  due  to  che- 
mical reactions,  depassivation  of  steel  due  to  critical 
chloride  ion  concentration,  cathodic  and  anodic  poten- 
tial due  to  concentration  polymerization  of  electrodes, 
and  flow  of  electric  current  through  the  electrolyte  in 
the  concrete  pores.  Practical  applications  of  the  model  were 
described  in  the  second  part  of  Bazant 's  paper  (1979b). 

2 . 2 Fiber  Reinforced  Concrete 

The  basis  for  the  concepts  of  fiber  reinforced  con- 
crete was  provided  by  the  early  research  of  Romualdi  and 
Batson  (1963a),  and  Romualdi  and  Mandel  (1964)  on  the  effects 
of  the  closely  spaced  reinforcing  wires  on  the  tensile 
strength  of  concrete.  They  found  that  for  the  same  volume 
of  reinforcing  steel,  the  tensile  strength  of  concrete  is 
significantly  increased,  if  the  spacing  between  the  wires  is 
reduced  to  less  than  10.2  mm  (0.4  in.).  This  finding  can  be 
explained  from  Griffith's  fracture  theory  which  states  that 
the  tensile  strength  is  inversely  related  to  the  length  of 
the  critical  flaws,  with  an  inherent  assumption  that  the 


critical  flaws  are  of  the  same  dimension  as  wire  spacing. 
It  was  observed  that  the  results  of  the  closely  spaced 
continuous  wires  could  be  achieved  by  mixing  short  lengths 
of  fibers  directly  in  the  concrete. 

The  concrete  properties  have  been  improved  to  such  an 
extent  by  the  addition  of  fibers  that  many  applications  of 
it  are  foreseen  which  were  considered  to  be  impractical 
before.  The  increasing  attention  the  fiber  reinforced  con- 
crete has  been  receiving  is  reflected  by  more  than  2,000 
references  reported  by  Hoff,  Fontenot,  and  Tom  (1976)  and 
Hoff  (1977,1979).  ACI  Committee  544  (1973)  has  published  a 
State-of-the  art  report  on  fiber  reinforced  concrete.  The 
guidelines  for  measuring  the  properties  of  fiber  reinforced 
concrete  have  been  discussed  by  ACI  Committee  544  (1978).  In 
this  section  some  of  the  aspects  of  fiber  reinforced 
concrete  are  briefly  reviewed. 

2.2.1  Various  kinds  of  fibers 

Fiber  reinforced  concrete  is  a composite  material 
consisting  of  cement-based  matrix  and  an  ordered  or  randomly 
distributed  short  fibers.  Many  different  kinds  of  fibers 
have  been  developed  for  use  in  concrete.  Fiber  types  can  be 
divided  into  two  main  groups:  those  with  moduli  higher  than 
the  cement  matrix,  such  as  steel,  glass,  asbestos,  and  those 
with  moduli  lower  than  the  cement  matrix,  such  as  nylon, 
polypropylene,  and  kevlar.  Steel,  glass,  kevlar,  and  poly- 
propylene are  the  most  frequently  used  fibers.  Due  to  the 


health  hazards,  the  use  of  asbestos  fibers  has  been  dis- 
continued recently. 

Steel  fibers  are  usually  of  the  straight  type.  But  to 
improve  the  bond  special  types  like  hooked-end,  crimped,  and 
Duoform  have  been  developed.  Though  steel  fibers  have 
high  strengths  and  are  found  to  produce  stiff  composites, 
they  only  increase  the  first  crack  strength  marginally.  If 
the  composite  is  designed  to  fail  by  fiber  fracture,  they 
increase  the  ultimate  strength  considerably.  Some  of  the 

applications  of  these  fibers  have  been  as  overlays  for 
roads,  pavements,  bridge  decks,  airport  run  ways,  and  indus- 
trial flooring.  Due  to  their  susceptibility  to  corrosion, 
further  studies  are  needed  for  their  application  to  off- 
shore environments. 

Glass  fiber  reinforced  concrete  is  being  considered 
as  a substitute  for,  and  improvement  of  asbestos  fibers. 
The  major  problems  are  fiber  breakage  due  to  instability 
and  the  surface  degradation  by  the  high  alkalinity  of  the 
hydrated  cement  paste.  While  alkali  resistant  (AR)  fibers 
are  more  stable  than  ordinary  E-glass  fibers,  the  degree  and 
severity  of  time-dependent  strength  degradation  and  their 
performance  in  wet  conditions  are  still  subject  to  specu- 
lation. One  of  the  important  property  improvements  obtained 
by  glass  fiber  inclusion  is  the  spectacular  increase  in  the 
impact  strength  compared  to  asbestos  fibers.  The  use  of 
high  alumina  cement  and  pulverized  fuel  ash  in  the  matrix 
can  improve  the  durability  of  these  fibers. 
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Both  polypropylene  and  kevlar  fibers  belong  to 
the  class  of  plastic  fibers.  They  have  high  tensile 
strengths,  but  their  low-moduli  precludes  reinforcing 
effects.  However,  their  high  elongation  capacity  (15-25 
percent),  makes  them  absorb  large  amounts  of  energy, 

enabling  the  composite  to  absorb  10  to  25  times  more  impact 
energy  than  the  plain  mortar  or  concrete.  The  most  exten- 
sive use  of  polypropylene  fibers  is  in  concrete  piles.  The 
superior  impact  resistance  properties  of  these  fibers  is  due 
to  a purely  mechanical  bond  between  the  fiber  and  concrete. 
Though  kevlar  fibers  have  not  been  used  in  the  concrete, 
their  future  applications  seem  to  be  promising  due  to  their 
similarity  to  polypropylene  fibers. 

2.2.2  Suitability  to  the  marine  environment 

Though  fiber  reinforced  concrete  has  many  superior 
structural  properties  than  plain  concrete,  its  suitability 
to  marine  applications  will  depend  on  the  stability  of  the 
fibers  in  a marine  environment.  Only  a few  of  the  many 
publications  on  fiber  reinforced  concrete  pertain  to  its 
marine  applications. 

Barab  and  Hansen  (1973)  described  the  performance  of 
steel  fiber  reinforced  breakwater  armor  units.  After  one 
winter  of  exposure  to  ocean  storms,  the  steel  fiber  units 
were  found  to  have  no  damage,  whereas  some  of  the  similar 
non-fiber  units  placed  at  the  same  time  sustained  some  crack- 
ing. Fiber  reinforced  units  showed  superior  impact 
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resistance  compared  to  the  ordinary  units.  Due  to  the 
limited  duration  of  marine  exposure,  no  comments  were  made 
regarding  durability. 

Based  on  their  tests  on  metal  fibers,  Rider  and  Hei- 
dersbach  (1980)  recommended  stainless  steel  fibers  for  marine 
applications.  Concrete  beam  specimens  with  three  different 
types  of  steel  fibers,  stainless,  carbon,  and  meltex 
(melt-spun  stainless)  were  exposed  to  circulating  sea  water  for 
up  to  1 2 months.  Flexural  strengths  were  determined  as  an 
indication  of  the  fiber  durability.  No  evidence  of  corro- 
sion-induced damage  was  evident  on  any  of  the  stainless 
steel  fibers  examined  by  electron  microscope,  whereas  the 
other  types  had  traces  of  corrosion. 

Reddy,  Gopal,  Nomikos,  and  Thiel  (1984)  studied  the 
strength  and  durability  of  marine  exposed  concrete  rein- 
forced with  three  different  kinds  of  fibers,  steel,  glass, 
and  kevlar  with  and  without  latex  modification.  Compress- 
ive, tensile,  flexural,  and  impact  strengths  were  measured 
at  various  periods  of  sea  water  exposure  as  an  indication  of 
structural  integrity.  Long  term  durability  was  studied  with 
alternate  exposure  to  sea  water  and  oven  drying.  Kevlar 
fibers  were  found  to  be  more  stable  than  the  steel  and  glass 
fibers.  The  addition  of  latex  was  found  to  significantly 
improve  the  durability  characteristics. 


2.2.3  Mechanics  of  fiber  action 


The  fibers  act  as  crack  arrestors  that  restrict  the 
growth  of  flaws  in  the  matrix,  controlling  them  from  en- 
larging under  stress  into  cracks  which  eventually  cause 
failure.  This  process  considerably  increases  the  static 
and  dynamic  properties  of  concrete,  its  toughness, 
ductility,  and  impact  resistance.  The  major  factors  which 
control  the  performance  of  the  composite  are  the  relative 
fiber-matrix  stiffness,  volume  content  and  other  physical 
factors  of  fibers,  and  the  fiber-matrix  interfacial  bond. 


Relative  stiffness 

For  an  efficient  stress  transfer  to  the  matrix,  the 
elastic  modulus  of  fibers  should  be  much  higher  than  that 
of  the  matrix.  Based  on  simple  volume  proportions,  the 
composite  cracking  strength  in  static  tension  was  obtained 
by  Hannant  (1976)  as: 
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M = modular  ratio  = — ^ber — 

matrix 

It  can  be  easily  seen  from  Eq.  2.2.1  that  only  high  modulus 
fioers  with  high  'M1  values  can  bring  about  an  improvement 
in  the  strength  of  the  composite.  Low  modulus  fibers  with 


' M ' less  than  unity  decreases  the  strength.  High  strength, 
high  modulus  fibers  impart  characteristics  of  strength  and 
stiffness  to  the  composite.  Low  modulus,  high  elongation 
fibers  are  capable  of  large  energy  absorption  characteristics, 
and  impart  a greater  degree  of  toughness  and  resistance  to 
impact  and  explosive  loading.  The  former  also  contributes  to 
these  dynamic  properties  but  to  a lesser  extent. 


Critical  volume  percent 

Based  on  the  properties  of  fibers  and  matrix  in  the 
composite,  a critical  volume  percentage,  vf(crit)/  can  be 
found  which  will  enable  the  fibers  to  carry  the  load 
the  composite  sustains  before  cracking.  By  satisfying  the 
force  equilibrium  before  and  after  the  fracture,  the  critical 
volume  percentage  can  be  expressed  as  follows  in  the  form 
given  by  Hannant  (1976): 
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ultimate  cracking  strength  of  matrix, 
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E ^ = elastic  modulus  of  fiber. 

If  the  volume  content  of  fibers  exceeds  the  V- , . , , the 

f ( crit) 

failure  is  usually  by  debonding  which  results  in  multiple 
cracking  of  concrete.  This  is  a desirable  situation  because 
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it  changes  a basically  brittle  material  with  a single 
fracture  surface  and  low  energy  requirement  to  fracture 
into  a pseudo-ductile  material  that  can  absorb  transient 
minor  overloads  and  shocks  with  little  visible  damage. 

Interfacial  bond 

The  fiber-matrix  interfacial  bond  plays  an  important 
role  in  the  post-fracture  behavior  of  composites.  High  bond 
strengths  are  required  for  improved  tensile  resistance  but  a 
lower  interfacial  bonding  would  improve  the  fracture  tough- 
ness and  impact  resistance  through  energy  dissipation  and 
damping  at  the  interfacial  discontinuities.  If  the  inter- 
facial bond  is  such  that  the  composite  failure  occurs  by 
pull-out,  then  the  matrix  becomes  the  principal  tensile  load 
carrying  element,  with  only  a modest  increase  in  tensile 
strength.  To  achieve  a truly  two-phase  composite  action, 
the  matrix  must  be  so  designed  as  to  transfer  load  to 
the  fibers  so  that  they  contribute  fully  to  the  composite 
strength.  The  fiber  physical  parameters,  length  and  dia- 
meter, can  be  so  chosen  that  the  failure  of  composite  would 
be  by  fiber  fracture,  when  the  full  bond  strength  is  deve- 
loped over  one  quarter  of  the  length,  i.e.,  j . For  this  to 
happen,  the  bond  force  developed  should  be  higher  than  the 
load  required  for  fiber  fracture,  i.e., 
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Hence,  a suitable  aspect  ratio  can  be  chosen  for  the  fibers 
depending  upon  the  requirements.  The  ultimate  composite 
strength  for  the  case  of  a failure  by  fiber  pull-out  is 


a = N t -r 
cu  d 


a = N a j. 
cu  fu 


(2.2.5) 

and,  for  the  case  of  fiber  fracture, 

(2.2.6) 

where  N is  the  number  of  fibers  intersecting  a unit  cross- 
sectional  area.  For  random  fiber  distribution,  'N1  is 
obtained  by  probabilistic  methods,  discussed  in  the  sub- 
section 2.2.4 


Other  factors 

The  above  discussion,  and  the  equations  are  only 
valid  for  the  most  ideal  case,  when  the  fibers  are  aligned 
in  the  load  direction  and  the  full  bond  strength  is  devel- 
oped as  found  by  pulling  tests  on  a single  fiber.  If  the 
fibers  are  aligned  randomly,  some  correction  factors  have  to 
be  applied,  which  will  considerably  reduce  the  theoretical 
strength.  Actual  fiber  bond  strengths  are  less  than  the 
ones  found  in  the  pulling  tests,  due  to  the  fiber-fiber 
interaction.  An  analogy  to  this  interpretation  is  the 
decrease  of  the  total  bearing  capacity  of  a group  of  foun- 
dation piles  when  the  piles  are  brought  closer  to  each 
other . 

Zollo  (1975)  described  extrusion  techniques  which 
enable  the  most  favorable  fiber  orientation  and  bond 


21 


characteristics.  Naaman  and  Shah  (1976)  described  the  pull-out 
mechanism  in  steel  fiber  reinforced  concrete.  They  observed 
that  the  peak  pull-out  load  of  a single  fiber  is  much  higher 
compared  to  that  of  a group  of  fibers.  The  mechanics  of 
fiber  reinforcement  has  been  discussed  by  Shah  and  Rangan 
(1971),  Swamy,  Mangat,  and  Rao  (1973),  Swamy  (1974),  Hannant 
(1978).  McKee  (1969)  has  presented  equations  for  the  fiber 
spacing  in  terms  of  the  fiber  diameter  and  the  volume 
content. 

2.2.4  Probabilistic  analysis  techniques 

As  the  Eqs.  2.5  and  2.6  indicate,  the  theoretical 
strength  of  the  random  fiber  composites  depends  on  the 
probabilistic  number  of  fibers  intersecting  a unit  cross- 
sectional  area.  Some  of  the  probabilistic  models  are 
discussed  in  this  section. 

Matsuishi  and  Iwata  (1982)  studied  the  probabilistic 
strength  of  steel  fiber  reinforced  concrete.  The  expected 
value  of  mean  strength,  with  and  without  considering  the 
fiber  pull-out  effects,  was  obtained  as  a product  of  the 
probability  of  fibers  intersecting  the  cross  section  and  the 
load  carried  by  the  fibers.  For  the  case  with  fiber  pull- 
out effects,  a separate  function  which  depends  on  the  fiber 
lengths  was  constructed.  The  probabilistic  number  of  fibers 
per  unit  area  of  the  composite  for  the  random  fiber  distri- 
bution is  as  follows: 
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The  major  limitation  of  the  study  was  the  non-inclusion  of 
fiber  interaction  effects  in  the  strength  calculations. 

Naaman,  Moavenzadeh,  and  McGarry  (1974)  assumed  a 
Poisson  distribution  for  the  statistical  distribution  of  the 
fibers  in  the  concrete  mass.  The  tensile  member  was  assumed 
to  be  made  up  of  a chain  of  a series  of  N consecutive  links, 
breaking  at  the  weakest  link.  It  was  assumed  that  the 
failure  mode  after  matrix  cracking  is  ductile.  The  analyti- 
cal model  predicts  the  casual  effects  of  given  dominant 
parameters,  e.g.  the  volume  fraction  and  the  aspect  ratio  of 
fibers,  the  bond  and  tensile  strength  of  matrix  on  the 
composite  strength,  derived  either  from  the  data  or  from  the 
assumptions  on  which  the  model  is  based.  The  fiber  inter- 
action effects  were  considered  in  the  model.  The  following 
relations  were  derived  for  the  number  of  fibers,  bond- 
strength,  and  the  composite  cracking  and  ultimate  strengths 
for  a one  link  chain  member: 
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where 


v cc  = average  composite  cracking  strength, 

acu  = average  composite  post-cracking  strength, 

a = a constant  determined  experimentally, 

A = cross-sectional  area  of  the  link, 

t = average  bond  strength  determined  from 
pull-out  test  on  a single  fiber, 


Eq.  2.2.9  for  the  fiber  interaction  effects  was  found  by 
curve  fitting  the  experimental  results  on  0.254  mm (0.01 
in.)-  diameter  fibers. 

2.2.5  Models  for  fiber  behavior 

Shah  and  Rangan  (1971)  considered  a two-phase  compos- 
ite materials  approach  to  develop  the  stress-strain  curve 
of  the  composite  incrementally  by  using  secant  modulus  and 
Poission  s ratio  of  the  two  components  at  different  strain 
levels,  i.e. , based  on  the  properties  of  concrete  and  fibers 
at  different  strain  levels,  the  composite  properties  are 
determined  at  the  corresponding  strain  levels.  The  secant 
modulus  for  the  random  fiber  concrete  is  written  as: 


and 


T 


reduced  bond  strength  due  to  fiber 
interaction  effects. 


E. 


RC 


(2.2.13) 


E, 


'9  0 


where 


secant  modulus  with  transversely  aligned 
fibers , 
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and 

E = secant  modulus  with  longitudinally 
aligned  fibers. 

The  expressions  used  for  EgQ  and  EQ  as  proposed  by  Tsai 
(1964)  and  Hansen  (1965)  are  rather  complicated.  The  stress 
is  given  by: 

aRC  = ERC  £C  (2.2.14) 

By  repeating  the  process  at  different  strain  levels,  the 

complete  stress-strain  curve  for  the  composite  is  developed. 

Tanigawa  and  Hatanaka  (1983)  developed  mathematical 
expressions  for  the  hysteretic-strain  behavior  of  steel 
fiber  reinforced  concrete  under  repeated  compression.  They 
observed  that  the  envelope  equation  developed  by  Popovics 
(1973)  for  plain  concrete  is  valid  even  for  the  steel  fiber 
reinforced  concrete.  Different  expressions  were  developed 
for  loading,  unloading,  and  reloading  portions  of  the 
hysteresis  curve  in  terms  of  volume  content  and  aspect  ratio 
of  fibers. 

Keer  (1981)  developed  general  equations  based  on  the 
ACK  (Aveston,  Cooper,  and  Kelly)  theory  (1971),  to  describe 
the  loading,  unloading,  and  reloading  behavior  of  a cracked 
fiber  cement  composite.  In  particular,  the  theoretical 
relationships  between  the  unloading  strain  and  the  resulting 
residual  strain  at  zero  load  and  subsequent  reloading 
modulus  were  examined.  This  formulation  would  be  useful  in 
analyzing  the  fiber  composites  subjected  to  a short  term 
overload  or  overstrain  sufficient  to  crack  it. 
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Stroeven  (1982)  developed  constitutive  relations  for 
a cracked  region  of  a composite  based  on  morphological 
observations.  The  stress  transfer  capability  of  a cracked 
region  in  plain  or  fiber  reinforced  concrete  was  analyzed  in 
terms  of  the  micromechanical  aspects  of  the  particle  inter- 
lock mechanism. 

2.2.6  Fracture  mechanics  methods 

The  concepts  of  fiber  reinforced  concrete  originated 
with  the  application  of  linear  elastic  fracture  methods  by 
Romualdi  and  Batson  (1963b)  to  find  that  the  strength  of 
concrete  in  tension  is  limited  by  internal  holes  and  micro- 
cracks. It  was  found  that  the  tensile  cracking  strength  of 
concrete  increases  in  proportion  to  the  inverse  square  root 
of  the  reinforcing  bar  spacing.  Since  that  time  fracture 
mechanics  has  been  applied  to  many  fiber  composite 
analyses . 

Hillerborg , Modeer,  and  Petersson  (1976)  developed  a 
method  in  which  fracture  mechanics  methods  are  introduced 
into  finite  element  analysis  by  means  of  a model  where 
stresses  are  assumed  to  act  across  a crack  as  long  as  it  is 
narrowly  opened.  The  crack  is  assumed  to  propagate  when  the 
stress  near  the  crack  tip  reaches  the  tensile  strength.  When 
the  crack  opens,  the  stresses  are  not  assumed  to  fall  to 
zero  at  once,  but  to  decrease  with  the  increasing  crack 
width.  The  analytical  and  experimental  strengths  of  both 
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reinforced  and  unreinforced  beams  were  found  to  compare 
well . 

The  suitability  of  the  J-integral  as  a fracture 
criterion  for  fiber  reinforced  concrete  has  been  studied  by 
Mindess,  Lawrence,  and  Kesler  (1977).  The  J-integral  is  a 
path-independent  energy  line  integral,  a one  parameter  average 
measure  of  the  elastic-plastic  field  near  a crack  tip.  It 
can  be  considered  as  the  potential  energy  difference  between 
two  identically  loaded  beams  having  differing  crack  lengths, 
which  may  be  interpreted  as  the  energy  available  for  crack 
extension.  By  testing  notched  beams  with  differing  amounts 
of  glass  and  steel  fibers  in  bending,  the  J-integral  para- 
meter KIC  (critical  stress  intensity  factor)  and  GIC 
(critical  strain  energy  release  rate)  were  evaluated.  JIC 
was  found  to  be  a superior  indicator  of  the  effectiveness  of 
fiber  addition  than  the  other  parameters. 

A similar  study  carried  out  by  Velazco,  Visalvanich, 
and  Shah  (1980)  concluded  that  the  R-curve,  a plot  of  stress 
intensity  factors  against  the  corresponding  crack  lengths, 
is  more  promising  than  the  other  methods.  Various  fracture 
mechanics  methods  studied  include  critical  stress  intensity 
factor,  J-integral,  critical  crack  opening  displacement, 
compliance  techniques  for  determining  the  slow  crack  growth, 
and  R-curve  analysis.  They  tried  to  identify  a fracture 
parameter,  independent  of  test  specimen  geometry,  which  can 
correctly  predict  the  effects  of  fiber  addition,  and 
concluded  that  the  R-curve  approach  is  the  most  promising. 
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This  approach  was  also  found  to  be  satisfactory  by 
Wecharatana  and  Shah  (1983). 

2.3  Constitutive  Models  for  Concrete 

This  section  focuses  on  various  past  and  present 
constitutive  models  used  in  the  analysis  of  concrete.  The 
newly  developed  endochronic  theory  is  evaluated  with  re- 
spect to  the  conventional  plasticity  theories.  Recent  advan- 
ces in  the  endochronic  theory  are  briefly  summarized. 

Due  to  the  lack  of  proper  experimental  data,  concrete 
structures  have  been  analyzed  in  the  past  on  the  basis  of 
linear  elastic-fracture  models.  Though  these  are  attractive 
in  terms  of  their  simplicity  and  required  computational 
efforts,  they  are  accurate  up  to  a small  fraction  of  the  fail- 
ure loads  and  result  in  conservative  designs.  Recent  experi- 
mental work  on  complete  stress-strain  curves  in  both  uniaxial 
and  more  complex  loading  conditions  has  enabled  the  develop- 
ment of  curve-fitting  and  nonlinear  elastic  methods.  Equiva- 
lent uniaxial  models  have  been  formulated  for  biaxial  and 
triaxial  loading  conditions.  Also,  perfect  and  work-harden- 
ing  plasticity  theories  were  developed  expressing  plastic 
strains  in  terms  of  assumed  yield  surfaces.  These  models 
predict  different  plastic  behavior  for  different  yield  sur- 
faces. To  overcome  the  ambiguity  associated  with  the  choice 
of  yield  surface  in  the  plastic  theories,  endochronic  i.e., 

viscoplastic  theories,  have  been  developed.  Excellent 
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summaries  of  various  methods  have  been  presented  by  Nilsson 
(1979),  Chen  and  Ting  (1980),  and  Chen  (1982).  These  methods 
can  be  catogerized  into  various  groups  as  follows: 

1)  curve-fitting  methods, 

2)  linear  and  nonlinear  elastic  theories, 

3)  perfect  and  work-hardening  theories, 

4)  visco-elastic  and  visco-plastic  theories, 
and 

5)  endochronic  theories. 

2.3.1  Curve-fitting  methods 

Based  on  the  analysis  of  experimental  data,  many 
stress-strain  relations  were  proposed  in  terms  of  the 
maximum  principal  stress  and  the  corresponding  strain. 
Popovics  (1970,  1973),  Darwin  and  Pecknold  (1977),  and  Wang, 
Shah,  and  Naaman  (1978a)  discussed  various  curve  fitting 
functions  such  as  polynomials,  exponential,  fractional  and 
statistical  functions.  Popovics  (1970)  discussed  a nume- 
rical approach  to  describe  the  complete  stress-strain  curve 
in  terms  of  the  maximum  compressive  stress  and  the  corres- 
ponding strain.  This  was  later  verified  by  the  same  inves- 
tigator (1973)  based  on  large  experimental  data. 

A rational  approach  was  proposed  by  Wang,  Shah,  and 
Naaman  (1978a).  Separate  functions  were  constructed  for 
the  ascending  and  descending  portions,  in  terms  of  four 
constants  which  were  determined  from  four  conditions  that  the 
curve  passes  through  the  origin,  the  peak,  a specified  point 
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after  the  peak,  and  the  curve  has  a specified  secant  modulus 
at  45  percent  of  the  compressive  strength.  A similar  method 
was  proposed  by  Sargin  (1971)  in  terms  of  a parabolic  curve. 

All  the  above  models  are  limited  to  uniaxial  cases 
only.  A widely  used  formula,  proposed  by  Liu,  Nilson,  and 
Slate  (1972),  valid  for  both  uniaxial  and  biaxial  cases,  has 
the  following  form: 


a 


1 + 


ae 

e 


+ 


(2.3.1  ) 


where 

a,£  = stress  and  strain  in  the  principal  stress 
direction, 

o ,e  = experimentally  determined  values  of  maximum 

p principal  stress  and  the  corresponding  strain, 

and 

a = initial  tangent  modulus  in  uniaxial  loading. 

For  the  case  of  biaxial  loading,  'a'  is  defined  as: 

a = - (2.3.2) 

(1  - va) 

in  which 

a = ratio  of  the  principal  stress  in  the  orthogonal 
direction  to  that  in  the  direction  considered, 

and 

v = Poisson's  ratio  in  uniaxial  loading. 

The  basic  concept  is  to  treat  the  biaxial  stress- 
strain  behavior  of  concrete  as  an  equivalent  uniaxial  rela- 
tionship. The  strain  increment  in  each  principal  direction 
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is  evaluated  solely  by  the  principal  stress  in  that  direc- 
tion and  the  corresponding  tangent  stiffness.  Poisson's 
ratio  is  assumed  to  be  constant.  This  model  has  the 
disadvantage  that  it  is  applicable  only  to  certain  simple 
structural  elements,  like  beams  and  panels,  and  cannot  model 
the  hydrostatic  pressure  sensitivity  or  the  volume  increase 
near  failure. 

2.3.2  Linear  elastic  models 

In  these  models,  linear  elastic  relations  in  terms  of 
E and  v are  used  to  model  both  uncracked  and  cracked  beha- 
vior. If  a principal  stress  exceeds  the  tensile  strength  of 
concrete,  a crack  is  assumed  to  form  and  the  concrete  is 
modeled  as  an  orthotropic  material.  Stress  in  the  normal 
direction  to  crack  is  redistributed  to  the  adjacent  members 
in  the  structure.  When  the  compressive  strain  exceeds  a 
prescribed  value,  concrete  is  assumed  to  be  crushed,  loosing 
its  strength  completely.  Some  applications  of  linear 
elastic  models  have  been  discussed  by  Ngo  and  Scordelis 
(1967),  and  Suidan  and  Schnobrich  (1973). 

2.3.3  Nonlinear  elastic  models 

In  linear  elastic  models,  the  elastic  moduli  G and  K 
are  assumed  constant.  But  this  is  unrealistic  in  view  of 
che  nonlinearity  of  concrete.  In  nonlinear  elastic  models, 
they  are  expressed  as  functions  of  stresses  or  strains, 
usually  in  terms  of  octahedral  components.  Depending  on  the 
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way  stresses  and  strains  are  related,  the  nonlinear  models 
are  classified  into  three  different  types:  Cauchy,  Hyper- 

elastic, and  Hypoelastic.  In  the  Cauchy  type  formulation, 
the  current  state  of  stress  is  expressed  as  a function  of 
current  state  of  deformation  as  follows: 

°ij  = Fij(Ekl)  (2.3.3) 

This  formulation  does  not  depend  on  the  history  of  stress  or 

strain,  and  may  lead  to  a positive  energy  creation. 

Hyperelastic  (Green)  type  are  based  on  the  assumption 
of  the  existence  of  a strain  energy-density  function  W,  such 
that 

aw 

ai j ~3e..  (2.3.4) 
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in  which  W is  a function  of  the  current  components  of  stress 
and  strain  tensors,  defined  as: 


This  ensures  that  no  energy  can  be  generated  through  load 
cycles,  and  laws  of  thermodynamics  are  always  satisfied. 

Hypoelastic  formulations,  also  called  incremental 
formulations,  are  used  to  describe  the  material  behavior  in 
which  the  state  of  stress  depends  on  the  current  state  of 
strain  as  well  as  on  the  stress  path  followed  to  reach  that 
state.  The  incremental  constitutive  relations  for  time- 
independent  materials  are  written  as: 


e . . 


W 


(2.3.5) 


(2.3.6) 
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in  which  <3  and  §kl  are  the  stress  and  strain  increment 
tensors,  and  are  the  current  stress  tensor  components. 

Recently,  incremental  stress-strain  relations  have 
been  formulated  for  a special  class  of  hypoelastic  materials 
in  which  the  function  ' F 1 is  assumed  to  depend  on  the  stress 
(or  strain)  tensor  itself.  Murray  (1979)  has  summarized 
these  models  in  a comprehensive  manner.  Gerstle  (1981a  and 
1981b)  presented  simple  formulations  for  biaxial  and  tri- 
axial  behavior  in  terms  of  variable  bulk  and  shear  moduli. 
Kupfer  and  Gerstle  (1973),  and  Darwin  and  Pecknold  (1977) 
presented  incremental  biaxial  stress-strain  relations  for 
concrete  in  terms  of  variable  tangent  moduli  and  Poisson 
ratio. 


2.3.4  Perfectly  plastic  models 

Under  multiaxial  stress  cases,  especially  under  tri- 
axial  compression,  concrete  was  observed  to  flow  like  a 


ductile  material  on  the  yield  or  failure  surface  before  reach- 
ing the  crushing  strain.  Plastic  or  irrecoverable  strains 
are  assumed  to  take  place  once  the  yield  surface  is  reached; 
they  are  obtained  from  the  normality  conditions.  The 
complete  stress-strain  relationships  are  developed  in  three 
parts:  1)  before  the  yield,  2)  during  plastic  flow,  and 

3)  after  fracture.  The  total  strain  is  expressed  as  a sum 
of  elastic  and  plastic  parts  as  follows: 

d£i]  = deij  * dEij 


(2.3.7) 
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The  elastic  strains  are  determined  from  Hooke's  law  or  any 
nonlinear  elastic  model 


d£ij  ■ Cijkl  dokl 


(2.3.8) 


in  which  C,  is  an  elastic  modulus  tensor,  and  da  is 

ijki  i j 

the  stress  increment  tensor.  The  plastic  strain  increment 
is  written  as  follows: 


= dx 


in  which  g(  a^  ) is 
usually  taken  as  the 
dated  flow  rule 


9g 

3a.  . 


(2.3.9) 
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called  the  plastic  potential  surface, 
yield  surface,  f.  This  gives  an  asso- 


dX 


3f_ 
3 a 


13 


(2.3.10) 


in  which  d*  is  a positive  scalar  factor  of  proportionality 
that  is  nonzero  only  when  plastic  deformations  occur. 

Some  of  the  yield  surfaces  considered  by  previous 
researchers  include:  1)  the  Von-Mises  criterion,  2)  Extended 

Von-Mises  criterion  (or  Drucker-Prager  criterion),  and 
3)  Coulomb  or  modified  Coulomb  criterion. 


2*3.5  Work-hardening  plasticity  models 

The  most  recent  of  the  plasticity  theories  are  the 
strain-hardening  or  work-hardening  formulations.  These 
theories  are  developed  to  account  for  the  experimental 
observation  that,  after  a limited  plastic  flow,  the  mate- 
rials harden  requiring  higher  stresses  to  cause  further 
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® training . A yield  surface,  termed  a loading  surface, 
combining  both  perfect  plasticity  and  strain-hardening,  is 
postulated  and  an  associated  flow  rule  is  used  for  the 
plastic  strains  before  fracture. 

According  to  this  approach,  the  stresses  under  oper- 
ating conditions  are  expected  to  be  in  the  initial  discon- 
tinuous range  such  that  the  concrete  behavior  can  be  charac- 
terized as  linear  elastic.  This  initial  discontinuous 
surface  is  the  limiting  surface  for  elastic  behavior.  When 
the  material  is  stressed  beyond  the  elastic  limit  surface,  a 
subsequent  new  discontinuous  surface  called  the  loading  sur- 
face is  developed,  which  replaces  the  initial  elastic 
surface.  If  the  material  is  unloaded  from  and  reloaded 
within  this  subsequent  loading  surface,  no  additional 
irrecoverable  deformation  will  occur  until  this  new  surface 
is  reached.  If  straining  is  continued  beyond  this  surface, 
further  discontinuities  and  additional  irrecoverable  defor- 
mations result.  Three  different  kinds  of  hardening  rules 
are  commonly  used:  isotropic,  kinematic,  and  mixed. 

Bazant  (1977)  discussed  the  numerical  stability  of  the 
plasticity  theories  used  in  finite  element  analysis.  Rud- 
nicki  and  Rice  (1975)  proposed  an  incremental  plastic  theory 
which  exhibits  inelastic  dilatancy,  sensitivity  to  hydro- 
static pressure,  and  internal  friction.  Plasticity  theories 
have  been  reviewed  in  considerable  detail  in  a recent  book 
by  Chen  (1982). 
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2.3.6  Viscoelastic  and  viscoplastic  models 

None  of  the  plasticity  based  formulations  can  model 
the  time-dependent  or  strain  rate-dependent  behavior. 

Various  viscoelastic  and  viscoplastic  models  have  been 
developed  to  take  these  two  effects  into  account.  Visco- 
elastic models  have  been  applied  to  model  the  creep  and 
age-dependent  behavior  of  concrete  by  Bazant  (1975). 

2.3.7  Endochronic  models 

In  plasticity  theories,  plastic  or  irrecoverable 
strains  are  assumed  to  form  after  the  yielding  takes  place. 
But  in  reality  they  may  occur  right  from  the  beginning  due 
to  internal  grain  rearrangement  like  in  soil  or 
concrete,  which  is  not  accounted  for  in  the  plasticity  theo- 
ries. Another  limitation  of  plasticity  theories  is  that 
there  is  a certain  ambiguity  associated  with  the  choice  of  a 
proper  yield  surface.  Loading  and  unloading  criteria, 
needed  to  determine  the  plastic  strains,  cause  unnecessary 
complications  in  the  numerical  computations. 

The  endochronic  theory  proposed  by  Valanis  (1971) 
overcomes  some  of  the  limitations  of  plasticity 
formulations.  The  theory,  derived  on  the  basis  of  irreversi- 
ble thermodynamics,  has  been  applied  to  predict  the  mechan- 
ical response  of  metals  under  complex  strain  histories. 
Jain  (1982),  provided  an  excellent  summary  of  the  thermodyna- 
mic derivation  of  the  endochronic  constitutive  equations. 
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The  basic  concept  underlying  the  endochronic  formu- 
lations is  the  characterization  of  the  material  strain  state 
in  terms  of  an  intrinsic  time,  which  is  a non— decreasing 
scalar  variable  that  depends  on  the  increments  of  strain  as 
well  as  time,  and  geometrically  represents  the  length  of  the 
path  traced  by  the  states  of  the  material  in  a strain-time 
space  of  suitable  metric,  i.e.,  the  memory  of  the  material 
is  expressed  in  terms  of  an  intrinsic  parameter.  The  early 
versions  of  the  theory  were  rate-independent;  time  was  not 
involved  in  the  intrinsic  time  terms.  The  distance  along  a 
strain  path  between  two  strain  states  P and  P'  in  Fig.  2.2, 
denoted  as  d?  , termed  the  intrinsic  time  measure,  is  defi- 
ned as  follows: 


2 


P . , de  . . de.  .. 

ljkl  13  kl 


(2.3.11) 


The  tensor  P^kl  is  a material  property  tensor  as  it  varies 
from  material  to  material.  The  stress  at  point  'P'  is  not 
determined  simply  by  the  strain  at  'P',  but  by  the  history 
of  the  strain  along  the  path  OP.  Valanis  termed  the  mate- 
rials for  which  the  stress  is  a function  of  the  strain  his- 
tory with  respect  to  an  intrinsic  time  scale,  as  endochronic 
(Greek:  endo=inner,  chronos=time ) , and  the  corresponding 
theory  as  endochronic  theory. 

To  account  for  the  experimental  observation  that  the 
materials  strain  harden,  the  intrinsic  time  scale  Z is 
defined  as  follows: 

dZ  = ; — >0;  Z *0  (2.3.12) 

f(S)  d? 
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Fig.  2.2  Path  in  Strain  Space 
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in  which  f(^)  is  a monotonically  increasing  material  prop- 
erty-dependent hardening  function.  The  intrinsic  time 
scale  Z is  representative  of  the  deformation  history  of  the 
material.  The  basic  stress-strain  relations  of  Valanis  are 

t 

presented  below: 

If  the  temperature  effects  are  neglected,  Eqs.  3.24 
and  3.25  of  Valanis  (1971)  can  be  simplified  to  a form: 

3_fi; 

3Z 


f 9eii 

r . . = 2 | y ( Z - Z * ) — dZ' 

3 { 


(2.3.13) 


Considering  a single  exponential  term  for  the  modulus, 


u(Z)  = y^e 


-aZ 


(2.3.14) 


then 


aij  = 2 yo  Ie_a(Z  Z'}  de^(z'>  ' 
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in  which  de^  = — i-L  dZ ' 
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(2.3.15) 

(2.3.16) 


Integrating  Eq.  2.3.15  following  the  Leibnitz  rule,  results 
in 


do  . 


de.  . 
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+ y2-  dZ  a.  • 

2Vn  13 


(2.3.17) 


‘0  “^0 

The  first  term  on  the  right  in  the  above  equation  resembles 
the  classical  elastic  term  and  the  second  term  that  of 
plastic  term.  The  form  of  the  above  equation  recognizes  the 
experimental  observation  that  plastic  strains  can  occur 
immediately  following  the  loading.  The  plastic  strain  term 
in  Eq.  2.3.17, 


, P 

d£i  j 
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(2.3.18) 
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is  similar  to  the  Prandtl-Reuss  relation  of  classical 
plasticity 

= (dx  )a±j  (2.3.19) 

Endochronic  theory  contains  these  relations  of  classical 
plasticity  as  a special  case.  But  this  formulation  is  con- 
tinuous and  ' dZ  1 , does  not  depend  on  any  postulations  or 
inequalities.  On  the  contrary,  dx  needs  the  postulation  of 
a yield  surface  and  the  loading  and  unloading  criteria. 

Recent  applications 

Endochronic  theory  was  first  applied  to  concrete  by 
Bazant  and  Bhat  (1976)  by  extending  the  formulation  for 
metals  developed  by  Valanis  (1971  ) to  include  many  nonlinear 
features  of  concrete  like  hydrostatic  pressure  sensitivity, 
inelastic  dilatancy  due  to  deviatoric  strains,  strain  hard- 
ening, and  softening.  The  model  was  able  to  predict 
stress-strain  diagrams  for  concrete  of  different  strengths, 
failure  envelopes  of  uniaxial,  biaxial,  and  triaxial  stress- 
strain  diagrams,  lateral  strains  and  volume  expansion, 
hysteresis  loops  for  repeated  compression,  strain  rate 
effects,  and  creep  behavior.  This  is  the  only  constitutive 
theory  which  can  model  all  the  above  nonlinear  features  of 
concrete  at  the  present  time.  Various  endochronic  constants 
were  determined  by  fitting  large  sets  of  experimental  data. 

It  was  indicated  that  those  constants  would  be  the  same  for 
any  normal  weight  concrete  with  compressive  strengths  in  the 
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range  of  3500  to  6000  psi.  The  hysteretic  behavior  of 
concrete  beams  was  predicted  by  Bazant  and  Bhat  (1977)  using 
the  above  model.  The  same  formulation  was  applied  to  pre- 
dict the  finite  element  behavior  of  cubes,  beams,  and  other 
concrete  structural  elements  by  Sorensen,  Arnesen,  and 
Bergan  (1977). 

Endochronic  theory  has  also  been  extensively  applied 
to  soils.  Bazant  and  Krizek  (1976)  proposed  an  endochronic 
model  for  the  liquefaction  of  sand.  This  formulation  was 
applied  by  Ansal,  Krizek,  and  Bazant  (1982)  to  predict 
seismic  response  of  earth  dams.  Wu  and  Wang  (1983)  used  an 
endochronic  theory  to  predict  the  behavior  of  sand  under 
static  loading  and  Wu  and  Sheu  (1983)  described  a method  to 
model  shear  hysteresis  behavior  of  sand.  Bazant,  Ansal,  and 
Krizek  (1982),  and  Valanis  and  Read  (1982)  discussed  various 
applications  of  endochronic  theory  to  soils.  Gopal  and 
Reddy  (1985),  and  Reddy  and  Gopal  (1985)  applied  endochronic 
theory  to  analyze  the  frozen  soil  behavior  to  a static 
loading  in  the  Arctic  offshore  region. 

Recent  improvements 

Although  capable  of  modeling  many  complex  phenomena, 
early  endochronic  formulations  were  unable  to  properly 
characterize  the  response  to  small  cyclic  stress 
oscillations  superimposed  on  large  static  stresses.  Those 
formulations  do  not  predict  a positive  energy  dissipation  or 
closed  hysteresis  loops,  and  yield  a reloading  path  that  is 
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less  steep  than  the  unloading  slope  for  small  load-unload 
oscillations.  Sandler  (1978)  and  Rivlin  (1981)  pointed  out 
that  these  theories  implied  a physically  impossible  condi- 
tion of  positive  energy  creation  violating  the  Drucker 
stability  postulate.  It  was  indicated  that  this  feature  may 
lead  to  unstable  solutions  in  wave  propagation  problems. 
Valanis  (1981  ) has  addressed  some  of  these  criticisms. 
Details  of  those  criticisms  are  not  discussed  here  as  both 
of  them  were  directed  on  the  early  formulations.  Instead, 
the  latter  formulations,  which  overcome  the  deficiencies  of 
earlier  models,  are  discussed. 

Plastic  strain  formulations 

Valanis  (1980)  redefined  intrinsic  time  in  terms  of 
plastic  strains,  instead  of  the  total  strain,  which  enables 
the  formation  of  closed  hysteresis  loops.  With  this  new 
definition,  dC  in  Eq.  2.3.11  is  rewritten  as 
2 

= Pijkl  dEij  dekl  (2.3.20) 

in  terms  of  the  plastic  strain  increments.  It  was  estab- 
lished for  the  one-dimensional  case  that  the  rate  of  dissi- 
pation at  the  onset  of  loading,  unloading,  and  reloading  is 
zero,  leading  to  closed  hysteresis  loops.  The  new  model 
implies  the  existence  of  a yield  surface  in  an  'a  priori* 
fashion.  Various  versions  of  classical  plasticity  theories 
were  shown  to  be  particular  cases  of  the  endochronic  formu- 
lation. This  model  was  applied  by  Valanis  and  Read  (1982)  to 
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describe  the  soil  response  to  cyclic  shear  loading,  hydro- 
static compression,  and  densif ication  under  cyclic  shear. 
Lin  and  Wu  (1983)  applied  the  new  endochronic  formulation  to 
obtain  the  cyclic  behavior  and  plastic  strain  wave  propa- 
gation in  titanium  alloys.  Some  of  the  recent  developments 
of  the  endochronic  theory  and  their  applications  were 
discussed  by  Valanis  and  Lee  (1982).  More  recently,  Valanis 
and  Lee  (1984)  proposed  an  endochronic  theory  of  cyclic  plas- 
ticity which  predicts  very  closely  the  experimental  behavior 
of  mild  steel  and  grade  60  steel  under  various  loadings. 


Jump-kinematic  hardening  formulation 

Bazant  (1978)  introduced  the  concept  of  'jump-kine- 
matic hardening1  in  endochronic  theory  to  obtain  closed  hys- 
teresis loops  that  ensure  positive  energy  dissipation.  The 
concept  of  a yield  surface  with  a moving  center  was  used. 
The  basic  definitions  of  inelastic  deviatoric  and  volumetric 
strain  increments  were  modified  to  include  the  hardening 
coefficients.  In  the  original  paper  by  Bazant  and  Bhat 
(1976)  they  were  defined  as  follows: 

s . . 

deV . = dz  (2.3.21  ) 

J 2G 


and 


de"  = dX  + — dZ* 
3K 


(2.3.22) 


where 
deVj  , de: 


dZ,  dZ' 


= deviatoric  and  volumetric  strain  increments, 
= intrinsic  time  increments. 
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dX  = dilatancy  parameter, 
s1j  = deviatoric  stress  tensor, 

and 


o = volumetric  stress. 

In  the  kinematic  hardening  model  Eqs.  2.3.21  and 
2.3.22  are  modified  as  follows: 


dZ 


dei'j  ■ 


2G 


and 


de " = dX  + ( a - a ) \p 


in  which 


dZ 

3K 


(2.3.23) 


(2.3.24) 


ct  . . , 

13 


and 


current  centers  of  the  deviatoric  and  volu- 
metric loading  surfaces  in  stress  space, 


P , i|)  = coefficients  (less  than  unity)  which  reduce 
the  rate  of  accumulation  of  inelastic 
strains  in  cases  of  unloading  and  reloading 

P and  ip  satisfy  a certain  inequality  so  as  to  assure  closed 

hysteresis  loops.  For  non-cyclic  loading  a—  and  a can  be 

set  to  zero,  and  P and  to  unity.  For  the  case  of  cyclic 

loading,  centers  of  the  loading  surface  must  jump  into  the 

extreme  stress  points  whenever  loading  changes  to  unloading, 

or  unloading  to  reloading.  This  can  be  explained  as  follows 

in  reference  to  Fig.  2.3: 


When  the  loading  changes  to  unloading,  the  centers  of 
the  loading  surfaces,  and  cx  f are  moved  to  the  extreme 
stress  point  in  the  loading  path,  i.e.,  to  point  1.  When 
the  unloading  is  changed  to  reloading,  the  centers  are  moved 
to  the  extreme  stress  point  in  the  unloading  path,  i.e.. 
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Fig.  2.3  Closed  Hysteresis  Loops  in 

'Jump-Kinematic'  Hardening  Model 
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point  2.  Between  the  loading  reversals  they  are  kept  con- 
stant. This  process  was  found  to  give  closed  hysteresis 
loops  even  for  an  infinitesimal  load-unload  cases.  More 
details  on  this  formulation  can  be  obtained  from  Bazant, 

Ansal , and  Krizek  (1982),  and  Bazant,  Krizek,  and  Shieh 
(1983) . 

2.4  Flexural  Analysis  of  Beams 

2.4.1  Previous  experimental  investigations 

Based  on  their  experimental  investigation,  Shah  and 
Rangan  (1971),  stated  that  1.25  percent  by  volume  of  steel 
fibers  increases  the  toughness  of  concrete  (area  under  the 
load-deflection  diagram),  by  as  much  as  20  times  that  for  plain 
concrete.  Relative  toughness  was  found  to  increase  with  the 
aspect  ratio  (length  to  diameter  ratio)  of  fibers.  Snyder 
and  Lankard  (1972)  discussed  the  factors  affecting  the  flexu- 
ral strength  of  steel  fiber  reinforced  concrete.  Fiber  size 
limits  were  proposed  for  good  workability. 

Hughes  and  Fattuhi  (1977a)  discussed  the  effects  of 
the  addition  of  various  fibers,  including  fibrillated  and 
monofilament  polypropylene,  as  well  as  round,  straight.  Duo- 
form,  crimped  and  hooked-end  steel,  upon  the  flexural  pro- 
perties of  a basic  concrete  matrix  at  three  different  ages, 

30,  90,  and  190  days.  The  beams  were  tested  at  a constant 

central  deflection  of  0.029  mm  per  minute  under  third-point 
loading.  Duoform  steel  fibers  were  found  to  be  very  effec- 
tive due  to  their  superior  bond  characteristics.  This  was 
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in  contradistinction  to  crimped  steel  fibers  that  reduce  the 
strength  and  toughness  if  present  in  large  quantities,  due 
to  the  shape-induced  stress  concentrations. 

The  present  techniques  of  fibrous  concrete  flexural 
testing  have  been  reviewed  by  Zollo  (1980).  Based  on 
extensive  testing  with  varying  shear  span  to  depth  ratios, 
fiber  and  matrix  types,  it  was  concluded  that  testing  of  4x4x 
12  in.  beams  would  be  adequate  to  determine  the  flexural  char- 
acteristics. Based  on  his  experimental  data,  and  compar- 
able data  from  other  sources,  Johnston  (1982)  proposed  a 
relationship  which  permits  the  flexural  strength  in  any 
testing  arrangement  to  be  predicted  from  the  values  obtained 
in  any  other  arrangement.  A relationship  was  developed  to 
predict  the  flexural  strength  in  any  testing  arrangement, 
based  on  the  strength  of  6x6x18  in.  specimens  tested  under 
third  point  loading. 

2.4.2  Previous  analytical  investigations 

Ngo  and  Scordelis  (1967)  were  the  first  to  apply 
finite  element  techniques  to  the  analysis  of  concrete  struc- 
tures. They  analyzed  simply  supported  singly  reinforced 
concrete  beams  for  artificially  introduced  crack  patterns 
using  constant  strain  triangular  elements.  Concrete  was 
modeled  as  a linear  elastic  material,  and  the  bond  between 
the  steel  and  concrete  was  represented  by  bond-linkage 


elements . 
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Park,  Kent,  and  Sampson  (1972)  analyzed  the  response 
of  reinforced  concrete  members  to  cyclic  loading.  A layered 
beam  element  approach  was  used  to  obtain  the  load-deflection 
behavior  of  beams  by  an  incremental  method.  For  a given 
curvature  increment  at  the  mid-span  section,  the  internal 
moment  developed  is  found  by  an  iterative  procedure.  This 
moment  is  prorated  to  different  sections  along  the  length  of 
the  beam  and  the  corresponding  curvatures  are  determined. 
Once  the  curvatures  at  all  sections  are  found,  the  central 
deflection  is  obtained  by  integrating  the  curvature  diagram. 
The  same  procedure  was  applied  to  the  analysis  of  steel 
fiber  reinforced  concrete  beams  by  Reddy  and  Gopal  (1984). 

Bazant  and  Bhat  (1977)  predicted  the  hysteretic  beha- 
vior of  reinforced  concrete  beams  using  the  endochronic 
stress-strain  formulation  developed  by  them  in  1976.  They 
used  the  layered  beam  element  approach  as  discussed  before 
to  formulate  the  section-equilibrium  equations.  Stiffness 
matrices  of  concrete  and  steel  were  assembled  as  discussed 
by  Suidan  and  Schnobrich  (1973).  A perfect  bond  was  assumed 
between  steel  and  concrete.  Finite  element  applications  to 
concrete  structures  analysis  have  been  discussed  by  Phillips 
and  Zienkiewicz  (1976),  Valliappan  and  Doolan  (1972),  and 
Yuzugullu  and  Schnobrich  (1973). 

Raj  agopalan,  Parameswaran,  and  Ramaswamy  ( 1 974)  pro- 
posed probabilistic  equations  for  the  ultimate  moment  carry- 
ing capacity  of  beams.  Probabilistic  coefficients  were 
determined  for  various  fiber  alignments.  Assuming  a uniform 
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tensile  stress  distribution  below  the  neutral  axis,  Hughes 
and  Fattuhi  (1977b)  proposed  equations  for  the  flexural 
strength  of  cement  based  beams  in  terms  of  the  aspect  ratio, 
pull-out  load  for  a single  fiber,  and  the  critical  fiber 
length.  These  equations  were  found  to  compare  well  with  the 
experimental  strengths  of  steel  and  polypropylene  fiber- 
reinforced  concrete  beams. 

Parimi  and  Rao  (1973)  presented  a method  to  calculate 
the  fracture  toughness  of  fiber-reinforced  concrete,  based 
on  a realistic  estimation  of  the  spacing  between  randomly 
placed  fibers.  The  analytical  expressions,  derived  by 
energy  methods,  compared  well  with  the  experimental  data  on 
cement  mortar  and  concrete  specimens  with  steel  fibers. 
Swamy  and  Mangat  (1974)  presented  a combined  crack  arrest- 
composite  materials  theory  to  predict  the  first  crack  and 
ultimate  flexural  strength  of  steel  fiber  reinforced 
concrete  beams.  They  assumed  the  fracture  process  to 
consist  of  two  stages:  initial  stage,  during  which  fibers 

debond  with  a slow  crack  propagation,  and  the  final  stage 
during  which  the  failure  occurs  due  to  the  unstable  crack 
propagation  when  the  fibers  pull  out  and  the  interfacial 
shear  stress  reaches  the  ultimate  bond  strength.  Wecha- 
ratana  and  Shah  (1983)  proposed  nonlinear  relations  between 
the  crack  opening  pressure  and  its  opening.  They  divided 
the  effective  crack  length  into  three  zones:  1 ) a traction 

free  zone,  2)  a fiber  bridging  zone,  and  3)  a matrix  process 
zone  resulting  from  aggregate  interlock  and  microcracking. 
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The  crack  closing  pressure  in  the  matrix  process  zone  was 
assumed  to  be  negligible,  and  ignored.  The  traction  free 
zone  is  the  one  in  which  fibers  are  pulled  out  completely. 
Hence,  only  the  fiber  bridging  zone  offers  resistance  to 
the  crack  propagation.  Analytical  results  from  linear  elas- 
tic fracture  methods  were  found  to  compare  well  with  the 
experimental  data  on  notched  beams. 

Additional  load-carrying  capacity  due  to  fibers 

Additional  moment  carrying  capacity  of  concrete  beams 
due  to  steel  fibers  in  the  presence  of  the  conventional 
steel  reinforcement  was  discussed  by  Henager  and  Doherty 
(1976).  The  strength  due  to  the  fibers  is  directly  added  to 
that  of  the  reinforcing  bars  to  obtain  the  theoretical 
ultimate  moment. 

The  recent  trend  in  fiber  reinforced  concrete 
construction  is  to  optimize  the  use  of  fibers  to  minimize 
the  costs.  Nagaraj  and  Dwarakanath  (1984)  proposed  a method 
to  estimate  the  capacity  of  partially  fibrous  concrete 
beams,  in  which  the  fibers  are  only  provided  in  the  tensile 
zone  where  they  are  active  and  no  fibers  are  provided  in 
compression  zone  where  they  are  inactive.  Their  method  is 
briefly  discussed  below. 

The  following  assumptions  were  made: 

1 ) the  fibers  are  assumed  to  be  uniformly  distributed  and 


randomly  oriented. 
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2)  flexural  tensile  stress  a , at  which  the  ductile 

stress  block  builds  up,  is  taken  equal  to  the  modulus 

of  rupture  of  the  plain  concrete  a , i.e.,  a,  = a 

m t m 

3)  strain  in  the  composite  corresponding  to  = a is 

t m 

considered  identical  to  that  of  the  matrix  corres- 
ponding to  its  modulus  of  rupture 

4)  strain  distribution  across  the  depth  of  the  section  is 
linear 

Fig.  2.4  shows  the  assumed  stress  and  strain  distribution 
across  the  depth  of  the  section. 

From  the  strain  distribution  diagram. 


'mut 
: cut 


i.e. , 
where 


1 

t 

n = 


D - nD  - kD 


= 1 - 


n 


D - kD 

(1  - k)  (t  - 1) 
t 


1 - k 


(2.4.1  ) 
(2.4.2) 


t = 


’cut 


'mut 


'cut: 

'mut 


= ultimate  strain  of  the  matrix, 

= ultimate  strain  of  the  composite, 

= tensile  strain  enhancement  factor. 


and 


k = neutral  axis  depth  factor, 


n = height  of  the  fibrous  tensile  zone  to 
achieve  a balanced  failure. 

Balanced  failure  is  the  failure  at  which  the  ultimate  compo- 
site strain  ecut  , and  the  ultimate  matrix  strain,  £mut 
will  be  reached  simultaneously  at  the  tensile  extremity  and 
the  unreinforced  concrete  at  the  interface  respectively. 
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Fig.  2.4  Partially  Fibrous  Beam  Section 
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From  the  stress  distribution  diagram,  the  equilibrium  of  the 
section  can  be  written  as  follows: 

o a 

amnD  + — (1-k-  n)D  = _£  k D (2.4.3) 

2 2 

By  solving  Eq.  2.4.3,  k can  be  written  as 


k = 


in  which 


1 + n 
1 + r 


(2.4.4) 


r = 


On  substitution  of  'n'  from  Eq.  2.4.2,  k becomes 


k = 


2 ~ T 
2 " E+  r 


(2.4.5) 


The  ultimate  moment  capacity  can  be  computed  as 


M = T1  [ n j + (1  - k-  n)  D + jkD  ] 

+T2  [ | (1  -k-n)  D+jkD]  (2.4.6) 

where 

T1  = CTm  B n D (2.4.7) 

and 

T2  = I amB  (1  ~ k - n)  D (2.4.8) 

Substituting  expressions  for  Tx  and  T2  in  Eq.  2.4.6,  and 
neglecting  the  second  order  terms  in  'n',  M can  be  written 
as : 


M 


6 a 


m 


(2.4.9) 


= 2 (1  + n - k), termed  the  moment  capacity 

enhancement  factor 


in  which 


B 
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The  tensile  strain  enhancement  factor,  ' t'  is  deter- 
mined from  a direct  tension  test  on  fiber  concrete  speci- 
mens. The  analytical  formulation  was  found  to  compare  very 
well  with  the  experimental  results  on  steel  fiber  reinforced 
concrete  beams  with  varying  fiber  contents. 


CHAPTER  III 

ENDOCHRONIC  STRESS -STRAIN  MODEL  FOR  CONCRETE 
3 . 1 Introduction 

A complete  stress-strain  curve  consists  of  an  ascen- 
ding part  during  which  the  stress  in  the  specimen  increases 
and  a descending  part  during  which  the  stress  decreases;  in 
both  regions  the  strain  continues  to  increase.  A complete 
stress-strain  curve  is  needed  to  predict  nonlinear  behavior 
accurately  for  extreme  loading  conditions,  e.g.,  impact  and 
earthquake  loadings.  The  improvement  in  the  concrete  per- 
formance due  to  the  fiber  addition  is  reflected  by  the 
larger  area  under  the  stress-strain  curve.  After  the 
concrete  cracks,  the  fibers  bridging  the  cracks  provide 
resistance  to  the  further  crack  growth  and  separation  of 
concrete  into  fragments.  Fiber  reinforced  concrete,  there- 
fore, behaves  like  a pseudo-ductile  material,  with  higher 
stresses  after  the  peak  in  unloading  portion  than  the  plain 
concrete. 

When  a concrete  specimen  is  loaded  in  a conventional 
testing  machine,  both  the  specimen  and  the  machine  deform. 

When  the  stress  in  the  specimen  reaches  a peak  value,  there  is 
a reduction  in  the  load  which  results  in  a sudden  release  of 
stored  strain  energy  by  the  machine.  This  released  strain 
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energy  causes  uncontrolled  failure  of  the  specimen.  As  a 
result,  a stable  descending  part  cannot  be  observed.  In 
order  to  avoid  uncontrolled  failure,  the  load  on  the  speci- 
men must  be  reduced  at  a sufficiently  fast  rate  to  control 
the  rate  of  strain  energy  release  of  the  testing  system. 

3 . 2 Techniques  for  Obtaining  Complete 
Stress-Strain  Curves 

There  are  basically  two  methods  of  obtaining  a com- 
plete stress-strain  curve.  They  are  as  follows: 

1 ) use  of  a closed-loop  testing  system,  and 

2)  in  parallel-testing  with  a coaxial  special  alloy 
steel  tube. 

3.2.1  Closed  loop  testing  system 

In  a closed  loop  testing  system,  any  experimental 
output  can  be  chosen  and  controlled  as  an  independent  varia- 
ble. The  experimental  parameter  is  matched  against  its  pro- 
grammed counterpart.  If  a difference  exists,  an  error 
signal  is  generated.  This  signal  modifies  the  movement  of 
the  loading  platen  (cross  head)  according  to  the  polarity  and 
the  magnitude  of  the  signal  so  as  to  minimize  it.  For  exam- 
ple, in  a uniaxial  compression  test,  the  circumferential 
strain  rate  can  be  kept  constant  to  obtain  a stable  descen- 
ding  part.  The  speed  at  which  the  testing  system  responds 
is  an  important  factor  in  obtaining  a stable  descending  por- 
tion of  the  curve.  It  depends  on  the  frequency  response  or 
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the  response  time  (electronic  as  well  as  hydraulic)  of  the 
system. 

Hughes  and  Fattuhi  (1977c)  used  this  method  to  obtain 
the  stress-strain  behavior  of  fiber  reinforced  concrete 
cubes.  A Losenhausen  EHRI  control  unit,  activated  by  three 
LVDT  tranducers,  was  used  to  maintain  a constant  rate  of 
deformation. 

3.2.2  Testing  with  a steel  tube  in  parallel 

In  this  method  concrete  cylinders  are  loaded  co- 
axially in  parallel  to  a steel  tube  whose  diameter  is  larger 
than  that  of  the  specimen.  The  tube  is  made  of  a heat  treat- 
ed special  alloy,  so  that  it  displays  a linear  elastic 
behavior  up  to  the  desired  maximum  strain.  It  is  so  designed 
that  the  load  on  the  system  keeps  on  increasing  even  when 
the  specimen  is  unloading.  Load  on  the  concrete  specimen  is 
determined  using  the  strain  and  elastic  modulus  of  steel 
tube.  Strain  in  the  concrete  is  assumed  to  be  same  as  that 
in  the  steel  tube. 

This  method  offers  an  advantage  over  the  other  method 
in  that  the  machine  characteristics  are  not  a determining 
factor  in  obtaining  a stable  descending  portion  of  the 
curve.  In  spite  of  its  simplicity  it  has  certain  limita- 
tions: the  size  of  the  concrete  cylinder  is  limited  to 
accommodate  the  additional  load  on  the  steel  tube,  and  the 
failure  of  the  specimen  cannot  be  observed. 
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Wang,  Shah,  and  Naaraan  (1978b)  used  this  method  to 
obtain  complete  stress-strain  curves  of  normal  and  light- 
weight concrete  in  compression.  Shah,  Gokuz,  and  Ansari 
(1981)  summarized  various  experimental  techniques  available 
for  obtaining  complete  stress-strain  curves. 


3.3  Stress-Strain  Relations 

The  endochronic  formulation  proposed  by  Bazant  and 
Shieh  (1978),  used  in  this  investigation,  belongs  to  the  family 
of  theories  with  jump-kinematic  hardening  criteria,  as  de- 
scribed in  Section  2.3.  As  the  details  of  this  formulation  are 
discussed  at  considerable  length  by  Bazant  and  Bhat  (1976)  and 
Bazant  and  Shieh  (1978),  only  the  salient  points  are  presented 
here. 

The  deformations  are  assumed  to  be  so  rapid  that  the 
material  is  modeled  as  a rate-insensitive  material  by  drop- 
ping the  time-dependent  intrinsic  terms.  The  elastic  prop- 
erties are  assumed  to  be  incrementally  isotropic  and  the 
stress-strain  relations  are  written  separately  in  terms  of 
deviatoric  and  volumetric  components  as  follows: 


deij  = 

deij 

+ deij 

n 

2G 

+ -ii  dz 

2G 

(3.1a) 

de  = 

deS 

+ deP 

= 22.  + 
3K 

dX  + ^iZ  ' + dX  ' 
3K 

(3.1b) 

where 

©i  j = e1j  ~ £ = deviatoric  components  of  the 

strain  tensor. 
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and 


s . . 

ID 


e 


a . . 

ID 


1 

Tkk  = 


volumetric  strain, 

Kronecker  delta, 

deviatoric  components  of  the 
stress  tensor  a 

volumetric  stress, 

bulk  and  shear  moduli, 

inelastic  dilatancy, 

shear  compaction. 


Z,  Z'  = intrinsic  times  for  distortion 
and  compaction  respectively. 

The  increments  of  intrinsic  times  are  defined  as  follows: 


(3.2) 

(3.3) 

d£  = — — — ; dn  = F(e,a)d£ 
f (n  ,e  ,cr) 

dc  ' = — £n!  ; dn  ' = H(a;ds' 

h(n’) 


and 


where 
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d? 
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1 
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22 
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The  increments  of  the  inelastic  dilatancy  and  shear 
compaction  are  defined  as  follows: 
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dX  = £(A  ) L ( A,e,c)d£  (3.4) 

and 

dA'  = £ * ( X * ) L*(  X ' r c f o ) d£  (3.5) 

where 


J2  = second  invariant  of  the  deviator  of 
the  tensor  which  follows  in 
parentheses , 

I1  = first  invariant  of  the  tensor  which 
follows  in  parentheses. 


f (n  ,£  ,cr  ) , F(e  ,cr ) 
h(n ' ) , H(a ) 
£(X),  L( A ,e ,a  ) 


and 


constants , 

distortion  hardening  and  softe- 
ning functions  respectively, 

compaction  hardening  and  softe- 
ning functions  respectively, 

dilatancy  hardening  and  softe- 
ning functions  respectively. 


£ / € 1 = non-decreasing  scalar  variables 
called  distortion  and  compac- 
tion measure  respectively. 


3.3.1  Strain-hardening  functions 

Inelastic  strain  increments,  whether  due  to  plastic 
slip  or  due  to  microcracking,  decrease  as  the  residual 
strain  increases.  This  is  due  to  the  fact  that  the  potential 
locations  of  stress  peaks  gradually  decrease  with  the  accumu- 
lation of  strain,  which  is  evident  more  in  the  cyclic  loading. 
As  c,  and  £ are  adopted  as  the  measures  of  the  accumulated 
inelastic  strain  (or  damage),  function  ' f ' and  'h'  must 
increase  as  £ and  t,'  increase.  By  analysis  of  test  data, 
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the  strain  hardening  functions  were  defined  as  follows  by 
Bazant  and  Shieh  (1978): 

6,n  + 87n2 

f(n*e,a)  = ( 1 + — ) F , (3.6) 

1 + F2/a? 

F = 1 + 

3 J2U)  (l+a9/n2) 

and 

h(n ' ) = 1 + 2_L  + ( III  )2.  (3.7) 

63  B4 

In  the  above  equations,  F3  is  a function  which  makes  the 
cyclic  stress  loops  to  contract  more  strongly  at  low  cyclic 
deviatoric  strains.  Function  F^  governs  the  strain- 
softening behavior  after  the  peak. 

3.3.2  Strain-softening  and  hydrostatic 
pressure  sensitivity 

If  the  functions  F(e,a)  and  H(a)  remain  constant, 
the  uniaxial  stress-strain  diagram  would  approach  an  asymp- 
tote of  positive  slope  in  which  case  the  failure  of  the 
material  will  not  be  modeled.  To  obtain  a gradual  decrease 
of  slope  on  approach  to  peak  stress,  the  inelastic  strain 
increments  must  be  increased,  obtained  by  an  increase  in  F 
and  H.  Inelastic  volume  change  due  to  hydrostatic  pressure 
is  taken  into  consideration  through  the  intrinsic  time 
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• r ^ 

variable,  ^ called  compaction  measure.  By  analysis  of 
test  data  the  following  softening  functions  have  been  found: 


F ( e , cr ) 
where 


F + f 
1 2 


(3.8) 


and 
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(J2U)  }4 


J23 


a22  + ^ J2  ^ * 


)3  , 


{ J9  ( e ) } 4 , 

f4  = ( " r > ' 

a + {J  (e)  r 


and 


Fj-  = a,,  a . (l+a,0a  . ) {- 

5 11  mm  12  ram 


J2(e)  * 

a13%in!!S  + {J2(  = > ^ 


where 

amin'amed'  and  Vax  are  the  medium, 

and  maximum  principal  stresses.  aQ  , a^  , ...  a22  are 
endochronic  constants. 


3.3.3  Dilatancy  and  compaction  functions 

The  dilatancy  functions  govern  the  microcracking  in 
the  concrete.  Since  microcracking  is  primarily  due  to  devia- 
toric  strains,  it  is  expressed  as  a function  of  J2 (e ) . The 
increments  of  the  dilatancy  parameter  dX , are  given  as  follows 
dX  = i (X  ) L ( X ,e  ,o  )dC  (3.9) 


where 


A (X  ) 


= 1 - 


and 


L ( X ,e  ,ct  ) = 


1 - 


C3  X 2 

£ / ( z + ( 

c1I1(a)  ^o  c22  + J2(e) 


C4 J2  ^ e ^ 2 

-=-= ) > 
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Since  the  inelastic  strains  are  due  to  the  formation  of 
microcracks,  the  incremental  elastic  moduli  should  decrease 
as  the  microcracking  increases.  The  shear  and  bulk  moduli  G 
and  K are  expressed  as  follows: 


and 


G = 


K = 


1 + Cg  X 2(1  + v ) 


1 + c5X  3(1  - 2v ) 


(3.10a) 


(3.1  Ob) 


In  the  above  equations  c^ , C2  ...  Cg  are  endochronic  con- 
stants and  v is  the  poisson  ratio. 

Inelastic  strains  due  to  hydrostatic  pressure  are 
modeled  by  a separate  intrinsic  term  called  shear  compac- 
tion. It  was  assumed  that  the  shear  compaction  will  not 
have  any  effect  on  the  moduli.  The  following  functions  are 
constructed: 


*'(*')  = 
L ' (A  ' ,e  ,cr  ) = 

and 


c6  < 1 - 4^  » 

1/3° 
qmin  ^3 

1 + (g3/c8)3 


0.93 


93  = ' c7  amin ' ~ 2 > 


in  which  Cg  , c^ , and  Cy  are  the  endochronic  constants. 


64 


3 . 4 Experimental  Investigation 
3.4.1  Design  of  concrete  mix 

The  design  of  a fiber  reinforced  concrete  mix  differs 
from  that  of  the  usual  concrete  in  that,  a special  conside- 
ration is  needed  for  the  amount  of  cement  in  the  mix  and  the 
size  of  the  aggregate.  ACI  Committee  544  (1973),  reported 
that  a water  cement  ratio  of  0.4  to  0.6,  and  a cement  con- 
tent of  3200  to  5530  N/c.  m.  (550  to  950  lb/c.  yd.)  are 
needed  to  ensure  a sufficient  quantity  of  cement  paste  to 
coat  the  fibers.  To  ensure  a uniform  dispersion  of  fibers, 
the  aggregate  size  should  not  exceed  19  mm  (0.75  in.)  and 
preferably  9.5  mm  (0.375  in.),  Hannant  (1978).  The  use  of 
larger  size  aggregate  may  result  in  balling  effects,  where 
the  fibers  form  clusters  or  balls  resulting  in  a poor 
consistency  of  the  mix.  A maximum  size  lime  rock  aggregate 
of  9.5  mm  (0.375  in.)  corresponding  to  No.  89  stone  was 
chosen  to  ensure  good  workability  and  uniform  fiber  distri- 
bution. Steel  fibers  of  diameter  0.51  mm  (0.02  in.),  length 
of  30.5  mm  (1.2  in.),  and  0.5  percent  by  volume  were  used  in 
the  mix.  The  fibers  were  of  the  type  "Dramix",  ZP  30/50, 
with  hooked  ends  collated  in  bundles  with  water  soluble 
glue,  manufactured  by  Bekaert  Steel  Wire  Corporation.  After 
the  addition  of  fibers  to  the  concrete  mix,  the  glue 
dissolves  separating  the  fibers  which  then  disperse  through- 
out the  mix.  The  following  mix  proportions  were  chosen, 
based  on  the  guidelines  of  Portland  Cement  Association  (PCA) 
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published  in  1979,  which  gives  a 28  day  compressive  strength 
of  31  MPa  (4500  psi). 

Table  3.1 

Mix  Proportions  of  Concrete 


Item 

Content 

water 

2270 

N/c. 

m. 

(390  lb/c.  yd.) 

cement 

4950 

N/c. 

m. 

( 850  lb/c.  yd. ) 

coarse  aggregate 

61  70 

N/c. 

m. 

( 1 060  lb/c.  yd. ) 

fine  aggregate 

8420 

N/c. 

m. 

(1447  lb/c.  yd.) 

super  plasticizer 

1 750 

cc/c 

.m. 

(59.5  fl.  oz./c.yd. 

steel  fibers 

0 . 5 percent 

by  volume 

3.4.2  Fabrication  of  specimens 

Fiber  concrete  specimens  of  diameter  76.2  mm  (3  in.) 
and  length  152.4  mm  (6  in.)  were  cast  in  cardboard  molds  in 
three  layers  with  external  vibration  and  internal  tamping 
with  a thin  rod.  The  rod  used  was  thin  enough  not  to 
disturb  the  fiber  orientation.  The  specimens  were  removed 
from  the  molds  after  24  hours,  and  cured  for  about  a week 
with  water  sprinkling.  They  were  then  subjected  to  circu- 
lating sea  water  exposure  of  varying  durations. 

3.4.3  Alloy  steel  for  coaxial  loading  jacket 

A heat  treated  special  alloy  tube  was  used  to  encase 
the  concrete  cylinders  for  the  stress-strain  testing.  The 
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jacket  is  of  SAE  4140  grade  steel,  heat-treated  to  a Brinnel 
hardness  of  445.  The  chemical  composition  of  the  alloy  is 
given  as  follows: 


Table  3.2 

Chemical  Composition  of  Alloy 


Item 

% volume 

Item 

% volume 

carbon 

0.38  - 0.43 

manganese 

0.75  - 1 .0 

phosphorus 

0.04 

sulphur 

0.04 

silicon 

0.20  - 0.25 

nickel 

0.01 

chromium 

0.80  -1.10 

molybdenum 

0.15  - 0.25 

3.4.4  Test  procedure 

Complete  stress-strain  curves  were  obtained  by 
testing  the  concrete  specimens  in  parallel  with  a coaxial 
steel  tube.  One  end  of  concrete  cylinder  was  capped  with 
sulphur  capping  material,  according  to  ASTM  C-617  standards. 
The  concrete  cylinder  was  then  placed  inside  the  steel  tube 
and  capped  with  hydrostone  at  the  other  end,  to  make  it 
flush  with  the  top  of  the  steel  tube.  Fig.  3.1  shows  the 
schematic  diagram  of  the  load  sharing  mechanism.  A foil 
type  strain  gage  of  350  ohms  resistance,  0.5  inch  length  was 
fixed  at  each  end  of  the  diameter  at  mid-section  of  the 
steel  tube.  The  strains  were  measured  using  a quarter 
bridge  circuit  in  a Vishay  digital  strain  indicator.  Model 
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Fig.  3.1  Load  Snaring  Mecnanism 
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The  tests  were  carried  out  in  a Forney  hydraulic 
testing  machine  of  capacity  1335  kN  (300  kips).  Fig.  3.2 
shows  the  experimental  set-up  with  the  steel  tube  in  the 
machine  and  the  strain  gage  connections.  The  rate  of 
straining  was  10  micro-strains  per  second,  which  corresponds 
to  approximately  30  psi  per  second  on  the  concrete  cylinder. 
This  rate  lies  within  the  loading  rates  of  20  - 50  psi  per 
second,  recommended  by  ASTM  C-39  standards  for  compressive 
strength  testing  of  concrete  cylinders. 

3.5  Determination  of  Endochronic  Constants 

Various  constants  involved  in  the  hardening,  soften- 
ing, and  dilatancy  functions  of  the  stress-strain  rela- 
tions were  found  by  curve-fitting  the  experimental  behavior. 
At  first,  the  values  of  the  constants  were  varied  manually  to 
obtain  a curve  as  close  to  the  experimental  one  as  possible. 
Then,  an  IMSL  (International  Mathematical  and  Statistical 
Library)  subroutine  "ZXSSQ",  was  used  to  minimize  the  sum  of 
the  squares  of  deviations  between  the  experimental  and  pre- 
' dieted  curve.  The  subroutine  gives  a local  minimum  of  the 
sum  of  squares  of  M real  functions  in  R real  variables  by 
using  a finite  difference  Levenberg-Marquardt  algorithm. 

The  ordinates  at  various  strain  values  were  used  as  the  real 
functions.  The  global  minimum  was  obtained  by  repeating  the 
optimization  process  with  the  new  constants  found  in  the 
previous  iteration. 
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Fig.  3.2  Experimental  Set-up  for  the 
Stress-Strain  Tests 
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3.6  Results  and  Discussion 

Figs.  3.3  and  3.4  show  the  stress-strain  curves 
obtained  for  the  plain  and  the  steel  fiber  reinforced  con- 
cretes at  various  periods  of  sea  water  exposure.  The  advan- 
tage of  testing  with  a steel  tube  is  evident  by  the  control- 
led crack  pattern  shown  in  Fig.  3.5,  in  contrast  to  the 
crushing  failure  when  tested  without  the  steel  tube,  Fig.  3.6. 
The  compressive  strength  of  concrete  is  almost  the  same  as 
that  found  in  a separate  series  of  compression  tests  indi- 
cating that  the  presence  of  steel  jacket  does  not  adversely 
affect  the  behavior  of  concrete  specimen.  Fig.  3.7  shows  the 
endochronic  curve  fit  obtained  for  the  steel  fiber  reinforced 
concrete.  The  constants  for  the  biaxial  and  triaxial  cases 
could  not  be  determined  due  to  the  lack  of  experimental  data 
in  those  load  cases.  Those  constants  were  assumed  to  be 
same  as  those  for  plain  concrete  reported  by  Bazant  and 
Shieh  (1978).  The  constants  obtained  are  as  follows: 


ao 

= 0.65, 

al 

= 0.68/f', 

c 

a2 

fc  *5 

- 1175  < 4650>' 

a3 

= 324000/ (f^)^ 

a4 

= 0.045, 

a5 

= 2160/ ( f ^ ^ 

a6 
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Fig.  3.3 


Stress-Strain  Behavior  of 
Plain  Concrete 


72 


Fig.  3.4  Stress-Strain  Behavior  of  Hooked-end 
Steel  Fiber  Reinforced  Concrete 
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t 


Fig.  3.5  Crack  Pattern  in  the  Failed  Specimens 
Tested  with  Steel  Tube 
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Fig.  3.6  Crack  Pattern  in  the  Failed  Specimen 
Tested  without  Steel  Tube 
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Fig.  3.7  Endochronic  Curve  Fit  for  Hooked-end 
Steel  Fiber  Reinforced  Concrete 
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rH 

nJ 

25.0, 

a15 

= 

1 .095, 

ai6 

= 

1 .21 6, 

ai7 

= 

0.055, 

ai8 

= 

0.94, 

a19 

= 

6300/ ( f ' ) 
c 

a20 

= 

14.0, 

a21 

= 

1000.0, 

a22 

= 

0.04, 

f 1 

a23 

= 

0.2  r , 

bl 

= 

9,1  * 7029 

b2 

- 

fc' 

C1 

= 

2.0/f^, 

C2 

= 

0.003, 

C3 

= 

0.50, 

C4 

= 

2.0, 

c5 

= 

250.0, 

c6 

= 

0.002, 

C7 

= 

1 . 052x1 0~6, 

C8 

= 

0.001  , 

ei 

= 

30.0, 

62 

= 

3500.0, 

P3 

= 

0.08, 

e4 

= 

0.23, 

Z1 

= 

0.001 5, 

Z2 

= 

0.0125, 

X 

o 

= 

0.007, 

= 

0.007, 

and 


v = 0.182. 

The  initial  slope  formula  used  was  that  given  by  ACI, 


E = 57,000  / f' 

o ' c 


(3.11  ) 


where 


= compressive  strength  of  concrete  (psi). 

Fig.  3.8  shows  the  predicted  cyclic  compressive 
stress-strain  behavior  under  uniaxial  loading.  The  curves 
were  obtained  by  a strain  control  procedure  using  the 
constants  determined  above.  The  predicted  closed  hyste- 
resis loops  prove  the  superiority  of  the  model  compared  to 
the  earlier  versions  as  indicated  by  Bazant  and  Shieh  (1978). 


Stress  (MPa) 
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Fig.  3.8  Predicted  Hysteresis  Behavior  Under 


CHAPTER  IV 

FLEXURAL  BEHAVIOR  OF  FIBER  REINFORCED 
CONCRETE  BEAMS 


4.1  Introduction 

The  load-deflection  behavior  of  steel  fiber  reinforced 
concrete  beams  is  predicted  based  on  the  endochronic  theory 
for  concrete  discussed  in  the  previous  chapter  and  a proba- 
bilistic model  for  the  fiber  action.  Analytical  predictions 
are  compared  with  the  experimental  data. 

4.2  Material  Models 

4.2.1  Endochronic  theory  for  concrete 

The  endochronic  formulation  described  in  Section  3.3  is 
used,  with  the  constants  for  steel  fiber  reinforced  concrete 
given  in  Section  3.6. 

4.2.2  Model  for  fiber  action 

Fibers  are  assumed  to  participate  only  in  the  tensile 
region.  Their  effect  in  compression  is  accounted  for  by 
curve-fitting  the  stress-strain  behavior  under  compressive 
loading.  A crack  is  assumed  to  form  normal  to  the  principal 
stress  when  it  exceeds  the  tensile  strength  of  concrete. 
Experimental  data  on  the  fiber  pull-out  force  and  the  crack 
opening  displacement  (COD)  for  steel  fibers  were  reported  by 
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Naaman  and  Shah  (1976),  and  Matsuishi  and  Iwata  (1982). 
Before  the  full  interfacial  bond  strength,  t , is  developed, 
the  fiber  force  is  linearly  proportional  to  the  crack  ope- 
ning and  then  it  gradually  decreases  with  a further  increase 
in  the  crack  opening.  When  the  bond  stress  exceeds  the 
limiting  value  t , fiber  force  is  found  from  the  following 
expression  proposed  by  Wecharatana  and  Shah  (1983): 


a = 

a ( 1 - - )2  (4.1  ) 

max  n 

max 

where 

a = 

fiber  stress, 

a = 

max 

fiber  stress  when  the  full  interfacial 
bond  strength  is  developed, 

Z 

T d 

n = 

crack  mouth  displacement. 

n = 

max 

maximum  pull-out  displacement  of  fiber, 
Z/4 

and 


Z ,d  = 

length  and  diameter  of  fiber  respectively. 

The  total  tensile  resistance  provided  by  the  fibers 
is  computed  by  considering  the  probabilistic  number  of 
fibers  which  intersect  a unit  cross-sectional  area  as 
reported  by  Matsuishi  and  Iwata  (1982)  and  Naaman,  Moavenzadeh, 
and  McGarry  (1970)  in  the  following  form: 


N = 
s 

where 

2vf 

2 (4.2) 

Trd 

ii 

number  of  fibers  intersecting  a unit 
cross-sectional  area, 

where 
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and 


Vf  = volume  content  of  fibers. 


4.2.3  Model  for  reinforcing  bars 

The  reinforcing  bars  are  assumed  to  behave  elasto- 
plastically  until  the  initiation  of  strain  hardening.  The 
strain  hardening  portion  is  modeled  by  the  following  method 
as  proposed  by  Burns  and  Siess  (1962): 


112 ( e — e , ) + 2 e -e  . f 

f = f { — + — — — C—  - 1.7)  > 

y 60(e  -£  , ) + 2 e -e  , f 
s sh  u sh  y 


(4.3) 


where 


f = stress  in  the  steel  corresponding  to  a strain  of  e 
s s , 

fy  = yield  stress  of  steel, 

strain  at  which  the  hardening  will  commence, 

and 


eu  = strain  at  ultimate  stress  of  steel. 

If  the  strain  exceeds  the  yield  value  in  one  direction  of 
loading,  the  stress-strain  curve  would  no  longer  be  elasto- 
plastic  due  to  the  Bauschinger  effect;  it  will  be  a function 
of  the  residual  strain  for  the  load  reversal.  Upon  load 
reversal,  the  following  expression  is  used  for  stress  between 
zero  and  yield 

-2.05  £ 0.129  £ 

fs  = fv[  1 - exp(  — ) + ] (4.4) 

sh  e sh 

where 

= effective  strain  hardening  strain 
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It  is  a function  of  the  load  history  and  the  initial 
plastic  strain  , defined  as  the  difference  in  strain 
between  two  sequential  points  of  zero  stress  immediately 
preceding  a given  load  cycle  expressed  as  follows: 


i _ 


sh 


sh 

1.38 


in- 


liE 


(4.5) 


When  the  stresses  are  greater  than  yield  stress  on  the 
reversal,  the  virgin  stress-strain  curve  is  used  by  substi- 
tuting for  e ^ in  Eq.  4.3.  A perfect  bond  is  assumed 

between  the  reinforcing  bars  and  concrete. 


4.3  Load-deflection  Relations  of  Beam 

A layered  beam  element  approach  is  used,  which  is 
similar  to  those  used  in  earlier  investigations,  e.g.  Park, 
Kent,  and  Sampson  (1972).  The  internal  moment  for  any  given 
curvature  increment  is  determined  by  satisfying  the  equilib- 
rium of  the  cross  section.  Each  section  is  divided  into 
layers  within  which  the  strain  is  assumed  to  be  constant. 

The  moment  developed  at  the  center  of  the  beam  for  any 
curvature  increment  is  prorated  to  different  sections  along 
the  length  to  determine  the  corresponding  curvatures.  The 
central  deflection  is  found  on  the  basis  of  the  curvatures 
at  different  sections.  The  description  in  this  section 
closely  follows  that  of  Bazant  and  Bhat  (1972),  except  for  the 
inclusion  of  fiber  effects  and  the  deflection  calculations. 
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4.3.1  Incremental  equilibrium  equations 

When  the  concrete  is  uncracked  it  is  modeled  as  an 
isotropic  material.  The  stress  and  strain  are  related  as: 

{ da  + da " } = [ D ] { de } (4.6) 

where 


da  , da"  = elastic  and  inelastic  stress  increments, 
D = material  property  matrix, 


and 


de  = strain  increment. 

For  uncracked  concrete,  D is  written  as  follows: 


“ “ 

D11  °12 

°13 

0 

0 

0 

• 

D22 

D23 

0 

0 

0 

[ D ] 

__ 

°33 

0 

0 

0 

(4. 

D.  - 
44 

0 

0 

D55 

0 

D66 

where 

D11 

= d22  = 

D33  = 

K + 

4 

3 G, 

D12 

= D13  = 

D23  = 

K - 

1 G' 

°44 

= D55  = 

D66  = 

2G, 

G,  K 

= shear 

and 

bulk 

moduli 

of 

concrete 

determined  from  the  endochronic  theory. 
Inelastic  stresses  are  determined  from  Eqs.  3.1a  and 
3.1b  as  follows: 

daVj  = 2G  de?j  + 3K  deP  (4.8) 

When  the  principal  stress  exceeds  the  tensile 


strength  of  concrete,  cracks  are  assumed  to  form  normal  to 
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the  maximum  principal  stress  direction.  After  cracking,  con- 
crete is  modeled  as  an  orthotropic  material,  with  zero  stress 
across  the  crack.  The  concrete  modulus  matrix  in  the  local 
coordinate  system  of  cracks  shown  in  Fig.  4.1  is  given 
as  follows: 


[ D* 


D22 


D23 

°33 


0 
0 
0 
a D 


44 


0 
0 
0 
0 
a D 


’55 


66 


where 

D22 

<N 

CM 

Q 

II 

D23 

= D23 

D33 

= °3  3 

and 


D 


12 


D11 

D13D21 

°il 


D 


31 


11 


/ 


(4.9) 


a = shear  transfer  factor. 

[ D ' ] matrix  is  obtained  by  a transformation  of  the  matrix 
[D]  to  account  for  the  loss  of  stiffness  normal  to  the 
crack  plane.  Shear  transfer  factor  is  needed  to  account  for 
the  frictional  force  along  the  crack  plane.  A value  of  0.5 
is  used  for  a , reported  to  be  a reasonable  value  by  Suidan 
and  Schnobrich  (1973).  The  cracked  concrete  modulus  matrix 
[ D ' ] is  transformed  back  to  the  global  coordinates  by  the 


Fig.  4.1  Local  Coordinate  System  for  Crack 
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matrix  [ L ] of  the  direction  cosines  of  the  principal 
stresses  as  follows: 

[ D"  ] = [ L ]T  [ D'  ] [ L ] (4.10) 

[ D"  ] is  used  in  Eq.  4.6  in  place  of  [ D ] when  a crack 
occurs . 

Fig.  4.2  shows  the  coordinate  system  for  the  beam. 

The  cross  sections  of  the  beam  are  assumed  to  remain  plane 

after  bending.  The  shear  strains,  e , e , and  £ 

xy  yz  zx  ' 

and  the  transverse  normal  strains,  e and  e , are  taken 

yy  zz 

into  consideration  in  the  analysis.  They  are  assumed  to  be 

constant  throughout  each  cross  section.  Both  curvature  and 

normal  shear  e , are  assumed  to  contribute  to  deflections, 
xz 

The  longitudinal  normal  strain  at  any  level  of  beam  z ^ 
can  be  written  as  follows: 

dexx.  = de  ~ (4.11  ) 

where 

dexx  = incremental  longitudinal  strain  at  a distance 
i z^  from  the  reference  axis, 

de  = incremental  longitudinal  normal  strain  at  the 
reference  axis, 

and 

dk  = incremental  curvature  of  the  beam. 

The  shear  strains,  and  e , and  the  corresponding 

shear  stresses  can  be  neglected  for  in-plane  loading  of  the 
beam  in  xz-plane,  Fig.  4.2.  Then  the  cross-sectional  equi- 
librium equations  can  be  written  in  the  following  form: 


dM 


Z z. 


l 


+ Z 


i 


(4.12) 
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Fig.  4. 


2 Coordinate  System  for  Beam 
- Bazant  and  Bhat  (1977) 
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dN 

X 

= Z 

Ax. 

da 

XX  . 

+ 

Z 

dFx. 

(4.13) 

i 

l 

i 

i 

l 

dNy 

= Z 
i 

da 

YYi 

+ 

Z 

i 

dF 

yi 

(4.14) 

dNz 

= Z 

V 

dozz. 

+ 

z 

dFz. 

(4.15) 

i 

1 

i 

i 

i 

dV 

= ' Z 

Ax. 

da 

xz  . 

(4.16) 

and 

dV 

11  1 

The  summation  is  carried  out  on  all  the  layers.  A 
'z 


X . 


V 


1 ^ 1 

and  Az  are  the  areas  of  the  layers  in  the  yz , xz,  and  xy 
i 

planes  respectively. 


dF_ 


dF. 


longitudinal  tensile  force  developed  by  the 
fibers  in  the  case  of  a crack  opening 


n Ns  7 Ef  dsxx. 

1 


tensile  force  developed  by  the  fibers  in  the 
y-direction  in  the  case  of  a crack  opening 


n Ns  I ^ Ef  deyy , 


and 


where 


tensile  force  developed  by  the  fibers  in  the 
z-direction  in  the  case  of  a crack  opening 


n Ns  ? d2  Ef  dczz. 

1 


n = bond  efficiency  factor  for  the  fibers, 
taken  as  0.8, 


and 

Ef  = Young's  modulus  of  the  fiber  material. 

By  substituting  the  expressions  for  da  , da  , da__  , 

xx  yy  Z Z 

and  daxz  in  terms  of  de,  dk,  deyy  , dezz  , and  dexz  into 
the  cross-sectional  equilibrium  equations,  4.12  - 4.16,  five 
simultaneous  force  equilibrium  equations  are  obtained, 
represented  in  a matrix  form  as  follows: 
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{dfC}  = [ RC  ] {dE}  - {df"}  , (4.17) 

c T 

{df  } = internal  force  vector  of  concrete 


{dM, 

dN 

, dN 

, dN 

, dV} 

X 

y 

z 

{de. 

dk. 

dE  , 

dE 

, dE 

yy 

zz 

' xz 

{df"}  = vector  of  inelastic  forces. 


and 

Q 

[ R ] = stiffness  matrix  due  to  concrete  alone- 

The  force  equilibrium  equations  for  the  reinforcing  steel 
can  also  be  written  as: 

{dfS}  = [ RS  ] {dE}  (4.18) 


The  total  forces  in  the  cross  section  are  obtained  as  a sum 
of  the  contributions  of  concrete  and  steel,  Eqs.  4.17  and 
4.18,  which  yields 

{df}  = [ R ] {ds } - {df"}  (4.19) 

where 

[ R ] = [ RC  ] + [ RS  ] 

Solution  techniques  for  Eqs.  4.19  have  been  discussed  at 
length  by  Bazant  and  Bhat  (1977). 


4.3.2  Determination  of  deflections 

The  bending  moment  at  the  center  of  the  beam  is  pro- 
rated to  two  other  sections,  and  the  corresponding  curvatures 
are  determined.  Then  the  central  deflection  is  found  by  the 
moment-area  method.  The  load  acting  on  the  beam  is  found 
from  the  central  moment.  The  computer  program  discussed  by 
Bazant  and  Bhat  (1977)  has  been  modified  to  include  the  fiber 


89 


effects  and  to  predict  the  central  deflections  from  the 
curvatures . 

4.4  Experimental  Investigation 

The  experimental  data  used  was  that  from  the  author's 
tests  on  fiber  reinforced  concrete  beams  exposed  to  circu- 
lating sea  water  and  that  reported  by  other  investigators. 

The  experimental  results  of  Hughes  and  Fattuhi  (1977a)  were 
used  to  check  the  computer  program. 

As  part  of  this  investigation,  steel  fiber  reinforced 
beams  made  of  the  same  mix  described  in  Section  3.4,  exposed 
to  circulating  sea  water  were  considered.  Details  of  the 
beams,  shown  in  Fig.  4.3,  are  dimensions  of  102x102x 
357  mm  (4x4x12  in.),  reinforced  on  both  the  compression  and 
tension  sides  with  steel  bars  of  diameter  6.4  mm  (#2  bars) 
and  provided  with  stirrups  of  diameter  3.2  mm  (#1  bars)  at  a 
spacing  of  45  mm  (1.8  in.).  They  were  tested  in  3-point 
loading  over  a span  of  305  mm  (12  in.)  at  a constant  deflec- 
tion rate  of  1.25  mm  (0.05  in.)  per  minute.  The  tests  were 
carried  out  in  a Tinius  Olsen  fixed  head  hydraulic  testing 
machine  of  capacity  267  kN  (60,000  lb.).  The  central  deflec- 
tion of  the  beam  was  measured  as  the  machine  head  movement 
to  an  accuracy  of  0.025  mm  (0.001  in.)  by  a mechanical 
deflection  gage  as  shown  in  Fig.  4.4.  Three  beams  of  each 
kind  were  tested  and  the  average  values  reported.  The 
tensile  strength  of  the  concrete  was  obtained  by  testing 
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Fig.  4.J  Reinforcement  Details  for  Test  Beams 
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Fig.  4.4  Experimental  Set-up  for  the 
Flexural  Tests 
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cylindrical  specimens  of  diameter  76.2  mm  (3  in.)  and  length 
152.4  mm  (6  in.)  under  split  cylinder  tension,  Fig.  4.5. 

The  experimental  load-deflection  curves  reported  by 
Hughes  and  Fattuhi  (1977a)  were  obtained  by  testing  beams  of 
dimensions  102x102x508  mm  (span  = 406  mm)  under  third  point 
loading.  The  specimens  were  reinforced  with  1 .5  precent 
by  volume  of  Duoform  (0.64x59  mm)  and  crimped  (0.51x49  mm) 
steel  fibers.  Stress-strain  curves  for  these  two  concretes, 
used  to  determine  the  endochronic  constants,  were  obtained 
by  testing  cubes  of  102  mm,  Hughes  and  Fattuhi  (1977c). 

4.5  Results  and  Discussion 

Figs.  4.6  and  4.7  show  the  experimental  load-deflec- 
tion (L-D)  behavior  of  plain  and  hooked-end  steel  fiber  rein- 
forced concrete  beams  respectively.  Ductility  imparted  due 
to  steel  fibers  is  evident  from  the  larger  area  under  L-D 
curve  for  steel  fiber  beams.  Average  pseudo-toughness  index 
values  for  plain  and  steel  fiber  concrete  beams,  4.5  and  7.0 
respectively,  very  well  reflect  on  the  improvement  due  to 
fiber  addition.  The  pseudo-toughness  index  is  defined  as  the 
ratio  of  area  under  L-D  curve  to  80  percent  of  peak  load 
on  unloading  branch  and  the  area  under  L-D  curve  to  first 
crack. 

The  comparison  between  the  experimental  and  predicted 
behavior  of  crimped,  Duoform,  and  hooked-end  steel  fiber 
concrete  beams  is  shown  in  Figs.  4.8,  4.9,  and  4.10 
respectively.  The  model  over-predicts  for  crimped  fibers. 
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Fig.  4 . b Experimental  Set-up  for  the 

Split  Cylinder  Tension  Tests 
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Fig.  4.6  Flexural  Behavior  of  Plain 
Concrete  Beams 
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Fig.  4.7  Flexural  Behavior  of  Hooked-end  Steel 
Fiber  Concrete  Beams 
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Central  Deflection (mm) 


Fig.  4.8  Predicted  and  Experimental  Flexural 
Behavior  of  Crimped  Steel  Fiber 
Concrete  Beams 
★ 

From  Hughes  and  Fattuhi  (1977a) 
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Central  Deflection (mm) 


Fig.  4.9  Predicted  and  Experimental  Flexural 
Behavior  of  Duoform  Steel  Fiber 
Concrete  Beams 

★ 

From  Hughes  and  Fattuhi  (1977a) 
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Fig.  4.10  Predicted  and  Experimental  Flexural 
Behavior  of  Hooked-end  Steel  Fiber 
Concrete  Beams 
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under-predicts  for  Duoform  fibers,  and  gives  almost  exact 
values  for  hooked-end  steel  fibers.  This  seems  to  be  due  to 
the  inadequate  modeling  of  fiber  action.  Duoform  fibers  are  a 
special  type  with  high  bond  strength  in  contradistinction  to 
crimped  fibers  in  which  the  stress  concentrations  cause 
matrix  cracking  and  disbonding  of  fibers.  The  bond  proper- 
ties assumed  in  the  model  are  more  relevant  to  the  hooked- 
end  steel  fibers  as  indicated  by  the  good  agreement  in 
Fig.  4.10.  In  the  ideal  case,  fiber  pull  out  load-  slip 
distance  curves  for  each  kind  should  be  incorporated  in  the 
model  which  would  lead  to  improved  predictions. 


CHAPTER  V 

ENDOCHRONIC  FINITE  ELEMENT  ANALYSIS  OF  A 
FROZEN  SOIL  FOUNDATION 


5 . 1 Introduction 

Most  of  the  offshore  drilling  concepts  in  Frontier 
areas  are  of  the  gravity  type,  e.g.,  artificial  islands  and 
concrete  gravity  platforms  in  relatively  shallow  water  and 
subject  to  large  ice  and  wave  forces.  Artificial  islands 
have  been  reported  to  be  more  economical  in  terms  of 
construction  and  maintenance  costs  in  water  depths  of  less 
than  20  m.  (66  ft.),  Robertson  (1983).  In  the  foundation 
design  of  these  structures,  which  are  close  to  the  shore, 
advantage  can  be  taken  of  the  high  strength  frozen-soil 
layers.  If  the  thermal  structure  is  not  disturbed  these 
layers  provide  a firm  base  for  foundation  over  a long  dura- 
tion of  time. 

A frozen-soil  is  a thermal  condition  of  the  ground  in 

which  the  soil  is  perennially  frozen.  It  occurs  in  the 

o o 

areas  where  the  soil  temperatures  of  0 C (32  F)  and  below 
persist  over  two  consecutive  winters  and  the  intervening 
summer.  It  consists  of  solid  mineral  particles,  ice,  liquid 
(unfrozen)  water,  and  gaseous  inclusions.  The  properties  of 
it  are  highly  dependent  on  the  volume  contents  of  various 
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constituents  and  temperature.  Blouin  et  al.  (1979)  reported 
the  occurrence  of  subsea  frozen-soil  layers  in  the  Prudhoe 
Bay,  Alaska. 


5.2  Stress-Strain  Relations 

Experimental  data  on  frozen  soils  obtained  by  many 
investigators  suggest  that  its  behavior  is  highly  dependent 
on  the  test  temperature  and  rate  of  loading,  Ladanyi  (1981), 
Parameswaran  (1980),  and  Parameswaran  and  Jones  (1981).  Since 
ice  is  the  main  internal  bonding  agent,  even  a slight  change 
in  the  temperature  could  change  the  degree  of  cementation, 
affecting  the  properties  in  an  adverse  manner.  Parameswaran 
(1980),  in  his  experiments  on  frozen-soil  samples,  found 
that  the  specimens  deform  elastically  in  the  early  stages  of 
deformation  when  strain  is  less  than  one  percent.  As  the 
strain  rate  increases,  the  specimens  show  an  upper  yield 
point,  and  fail  essentially  in  a brittle  manner.  The  fail- 
ure stress  increased  with  increasing  rates  of  strain  and 
decreasing  temperatures. 

To  model  the  strain  rate  effects  and  temperature, 
stress-strain  relations  for  concrete  Eqs.  3.1a  and  3.1b  are 
modified  in  a simplistic  manner  as  follows: 


de . . = 
13 
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(5.1b) 


The  thermal  effect  is  incorporated  into  the  incre- 
mental volumetric  strain.  This  enables  higher  strain  incre- 
ments with  higher  temperatures  for  the  same  volumetric 
stress  increment.  The  thermal  strains  in  frozen  soil  are 
regarded  as  inelastic  in  view  of  their  irreversible  nature. 
The  intrinsic  times  dZ  and  dZ  ' are  redefined  as  follows 
to  include  the  time  effects: 

(dZ)2  = 

and 

2 

(dZ')  = 

where 

dt  = 
e = 

T1  = 

fcl  ' *2  = 

0 

and 

X # tk  = 
a d 

This  formulation  is  only  valid  for  negative  temperatures. 
The  time  hardening  function  for  creep,  t1  , makes  the  rate 
of  creep  decrease  with  time,  which  is  commonly  observed  in 
the  experiments.  This  method  is  the  same  as  that  applied  to 
concrete  by  Bazant  and  Bhat  (1976). 
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5.3  Experimental  Data 

The  experimental  data  used  in  this  investigation  is 

that  obtained  by  Parameswaran  (1980),  and  Parameswaran  and 

Jones  (1981)  from  tests  on  cylindrical  specimens  of  diameter 

50.8  mm  (2  in.)  and  length  140  mm  (5.5  in.)  with  20  percent 

3 3 

moisture  content,  and  a bulk  density  of  20kN/m  (127  lb/ft  ). 
The  tests  were  carried  out  at  different  temperatures  and 
strain  rates.  The  uniaxial  stress-strain  curve  was  obtained 
by  testing  the  specimens  at  a temperature  of  -6°C  (21.2°F)  at 
a strain  rate  of  3.98x10  ^ per  second.  The  triaxial  curve 
was  obtained  at  a temperature  of  -1 0°C  (14°F)  at  a strain 
rate  of  7.7x10  ^ per  second  and  a hydrostatic  confining  pres- 
sure of  0.1  MPa  (15  psi.) 

5 . 4 Determination  of  Endochronic  Constants 

The  endochronic  constants  are  determined  as  described 
in  Section  3.5.  Figs.  5.1  and  5.2  show  the  curve  fits 
obtained  for  the  uniaxial  and  triaxial  stress-strain  curves. 
The  lack  of  a good  fit  for  the  triaxial  curve  shows  the 
limitation  of  the  model.  The  curve  fit  could  be  improved  by 
constructing  some  other  functions  to  take  into  account  the 
melting  of  intergranular  ice  at  high  pressures;  this  re- 
quires a stronger  dependency  on  the  temperature  and  hydro- 
static pressure.  Various  constants  obtained  are  as  follows: 


a 


o 


0.52 


0.41/f 


c 


252,000/ ( f ^ ) 


4 


750  ( 
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Fig.  5.1  Uniaxial  Stress-Strain  Curve  Fit 
for  Frozen  Soil 

★ 

From  Parameswaran  (1980) 
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Fig.  5.2  Triaxial  Stress-Strain  Curve  Fit 
for  Frozen  Soil 

* 

From  Parameswaran  and  Jones  (1981) 
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• 

5.5  Finite  Element  Analysis  of  Structural  Response 


A finite  element  computer  program  has  been  developed 
based  on  the  endochronic  stress-strain  behavior  described  in 
Section  5.2.  It  can  account  for  the  time-dependent  creep 
and  thermal  effects.  Constant  strain  axisymmetric  elements 
are  used,  but  the  two-dimensional  problems  can  also  be  anal- 
yzed by  deleting  the  e term  from  the  element  equations.  For 

9 

any  given  load  increment,  the  nodal  displacements  and  the 
element  stresses  are  found  by  an  iterative  procedure.  For 
any  constant  set  of  loads,  the  creep  response  can  be  found  by 
incrementing  the  time  'dt'  in  Eqs.  5.2a  and  5.2b.  The 
computer  program  developed  in  this  section  closely  follows 
that  of  Bazant  and  Shieh  (1978).  The  accuracy  of  the  program 
is  checked  with  the  results  presented  for  a strip  footing  by 
Haidar,  Reddy,  and  Arockiasamy  (1982). 

The  foundation  response  to  the  static  weight  of  a 
typical  Arctic  offshore  island  is  obtained  by  the  computer 
program  discussed  above.  Details  of  the  island  as  reported 
by  Robertson  (1983),  are  a diameter  of  107  m.  (350  ft), 
finished  elevation  of  4 m.  (13  ft)  above  the  mean  sea  level, 
water  depth  of  3.35  m.  (11  ft),  and  approximately  120,035  m3 
(157,000  yd  ) of  gravel.  Fig.  5.3  shows  the  finite  element 
discretization  of  the  foundation.  All  the  nodes  on  the  two 
sides  and  the  bottom  are  assumed  to  be  supported  on  rollers. 
The  soil  is  assumed  to  have  a constant  temperature  over  the 
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Fig.  5.3  Finite  Element  Discretization 
of  the  Foundation 
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entire  depth  of  the  foundation.  If  the  temperature  profile 
is  known,  exact  temperatures  can  be  specified  for  each  ele- 
ment separately. 

5.6  Results  and  Discussion 

Fig.  5.4  shows  the  deflection  profile  at  -6°C  (21.2°F) 
temperature  at  various  times.  Figs.  5.5  and  5.6  show  the 
deflection  profiles  at  different  temperatures  and  times. 

As  expected,  the  deflections  are  larger  at  higher 
temperatures.  The  strength  and  the  initial  slope  of  the 
stress-strain  curve  of  the  frozen  soil  decrease  with  an 
increase  in  the  temperature.  This  is  clearly  verified  by 
the  higher  initial  and  creep-induced  deflections  for  higher 
temperatures.  With  a decrease  in  the  temperature,  the 
properties  of  the  frozen  soil  reach  steady  state  values, 
which  will  not  change  considerably  even  with  further 
temperature  reductions.  This  is  well  reflected  by  the 
proximity  of  the  deflections  at  -6  C and  -1 5 C temperatures. 

The  extension  of  the  endochronic  constitutive 
modeling  to  include  creep  and  thermal  effects  seems  reason- 
able, but  needs  to  be  studied  in  greater  detail  with  actual 
field  data  and  compared  with  other  analytical  approaches. 
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Fig.  5.4  Deflection  Profiles  at  Different  Times 
at  Soil  Temperature  = -6°C 
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Fig.  5.5  Deflection  Profiles  at  Time  = 0 
for  Different  Soil  Temperatures 
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Fig.  5.6  Deflection  Profiles  at  Time  = 2 years 
for  Different  Soil  Temperatures 


CHAPTER  VI 

SUMMARY  & CONCLUSIONS 


This  work  has  been  an  investigation  on  the  applica- 
tion of  endochronic  theory  to  the  analysis  of  marine  exposed 
fiber  concrete  structural  components  in  conjunction  with  a 
probabilistic  method  to  account  for  the  fiber  effects  and  an 
extension  of  the  endochronic  formulation  by  including  the 
thermal  and  creep  effects.  This  chapter  summarizes  the  find- 
ings of  this  investigation  and  scope  for  further  research 
in  this  area. 

6 . 1 Summary 

Initially,  a brief  review  of  the  present  state  of  the 
art  is  made  on  various  topics  pertaining  to  this  investi- 
gation, i.e.,  marine  concrete  behavior,  fiber  reinforced 
concrete,  constitutive  models,  and  concrete  structural 
analysis  techniques.  It  was  observed  that,  although  consid- 
erable work  has  been  carried  out  on  developing  models  for  the 
fiber  concretes  in  cyclic  loading  and  the  methods  to  predict 
the  additional  load  carrying  capacity  due  to  fibers,  not  much 
has  been  done  in  applying  these  in  incremental  structural 
analysis  programmes.  Most  of  these  theories  tend  to  be 
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based  on  the  particular  experimental  data  of  the  investi- 
gator and  may  not  be  realistic  for  other  fiber  types. 

The  endochronic  formulation  for  concrete  seems  to  be 
the  most  comprehensive  of  all  the  constitutive  theories 
capable  of  modeling  the  experimentally  observed  nonlineari- 
ties of  concrete,  i.e.,  strain  hardening,  softening,  dila- 
tancy  due  to  internal  microcracking,  and  hydrostatic  pres- 
sure sensitivity.  The  only  variable  which  needs  to  be  speci- 
fied to  describe  the  behavior  of  any  normal  weight  plain 
concrete  is  its  compressive  strength.  It  was  reported  that 
the  predicted  concrete  stress-strain  behavior  based  on  its 
compressive  strength  compares  well  with  the  experimental 
data  and  with  that  obtained  by  the  expression  of  Popovics 
(1973).  The  formulation  for  concrete  has  excellently  pre- 
dicted the  hysteretic  flexural  behavior  of  beams,  Bazant  and 
Bhat  (1977).  These  two  aspects  of  the  theory  prove  its  pre- 
dictive capabilities  in  contrast  to  the  usually  misconstrued 
belief  that  it  is  a just  another  curve-fitting  model. 

The  recently  developed  endochronic  models  which  over- 
come the  deficiencies  of  the  earlier  formulations  are  eval- 
uated along  with  some  of  their  applications.  Bazant  (1978) 
introduced  the  concept  of  jump-kinematic  hardening  into  the 
endochronic  theory  to  achieve  closed  hysteresis  loops. 
Valanis  (1980)  achieved  the  same  by  redefining  the  intrinsic 
time  in  terms  of  plastic  strains  instead  of  the  total 
strain.  Internal  constraints  in  these  two  formulations 
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require  an  additional  computational  effort,  but  still  pre- 
serve the  continuous  modeling  feature  of  the  theory. 

Complete  stress-strain  curves  were  obtained  by  testing 
the  concrete  specimens  in  parallel  with  a coaxial  special 
alloy  steel  tube.  This  technique  was  found  to  produce  a stable 
descending  portion  of  the  curve  irrespective  of  the  machine 
characteristics.  By  comparing  the  compressive  strengths  with 
those  obtained  in  separate  compression  tests  without  steel 
tube,  it  was  concluded  that  the  concrete  behavior  is  unaffec- 
ted by  the  presence  of  steel  tube.  Similar  conclusion  was 
made  by  Wang,  Shah,  and  Naaman  (1978a)  by  comparing  the  ascen- 
ding portion  of  many  curves  obtained  with  and  without  steel 
tube.  Endochronic  material  constants  were  determined  by  curve- 
fitting the  experimental  data  using  an  optimization  process. 

The  cyclic  behavior  of  steel  fiber  concrete  under  uniaxial  com- 
pression was  predicted  based  on  the  constants  determined 
from  a data  without  cycling.  Although,  it  could  not  be  com- 
pared with  any  experimental  curve  due  to  nonavailability, 
the  predicted  cyclic  behavior  seems  to  be  realistic. 

The  flexural  behavior  of  steel  fiber  reinforced  con- 
crete was  predicted  based  on  the  endochronic  formulation  for 
concrete  together  with  a probabilistic  model  for  the  tensile 
contribution  of  fibers.  The  fiber  model  is  based  on  the 
experimental  data  reported  by  various  investigators  on 
straight  and  hooked-end  steel  fiber  concretes.  The  major  fac- 
tors involved  are  the  probabilistic  number  of  fibers  inter- 
secting a unit  cross-sectional  area,  limiting  stress  at 


which  fiber  debonding  occurs,  and  the  frictional  shear  stress 
after  debonding.  The  present  model  is  not  valid  for  the 
cyclic  loading  cases  since  the  tensile  zone  fiber  behavior 
in  unloading  and  reloading  is  not  considered.  The  analyti- 
cal behavior  was  compared  with  the  experimental  data  for 
three  different  fiber  kinds,  hooked-end,  crimped,  and  Duo- 
form.  As  explained  in  Section  4.5,  the  model  is  quite  accu- 
rate for  hooked-end  fibers,  and  for  the  others  it  differs  by 
as  much  as  10  to  15  percent.  This  is  due  to  the  reason  that 
the  experimental  data  on  those  two  fiber  types  are  not 
incorporated  into  the  model. 

The  endochronic  theory  presents  a viable  method  for 
incorporating  various  features  of  material  behavior  into  the 
constitutive  equations.  In  Chapter  5,  a simple  method  was 
proposed  to  include  the  thermal  and  creep  behavior  of  mate- 
rials like  ice  and  frozen  soil.  The  response  of  a typical 
Arctic  island  foundation  consisting  of  frozen  soil 
was  analyzed  by  a finite  element  analysis  program  developed 
with  the  above  constitutive  law.  Endochronic  constants  for 
the  frozen  soil  were  determined  by  curve-fitting  uniaxial 
and  triaxial  stress-strain  curves  at  different  temperatures 
and  strain  rates.  The  results  of  this  program  have  been 
found  to  compare  favorably  with  those  of  ADINA  (Automatic 
Dynamic  Incremental  Nonlinear  Analysis)  as  shown  by  Gopal 
and  Reddy  (1985). 

Appendix-A  outlines  the  various  numerical  schemes 
used  in  the  present  investigation. 
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6 . 2 Suggestions  for  Further  Studies 

The  following  suggestions  are  made  for  future 

research  efforts  in  this  area: 

1 ) A model  should  be  developed  for  random  fiber  concretes 
after  cracking  under  cyclic  loading  by  adapting  similar 
models  for  continuous  fibers.  This  kind  of  a model  is 
needed  to  study  the  behavior  of  fiber  reinforced  concrete 
structural  elements  subjected  to  dynamic  loading. 

2)  It  is  desirable  to  simplify  the  present  endochronic  formu- 
lation for  concrete  to  enable  easy  determination  of  the 
constants  for  new  sets  of  data. 

3)  In  view  of  its  capabilities,  endochronic  theory  provides 
a vast  scope  for  developing  accurate  fracture  analysis  of 
fiber  concretes.  Most  of  the  previous  studies  have  been 
based  on  linear  elastic  methods  which  are  not  realistic 
for  concrete  due  to  its  nonlinearities.  Ductility  due  to 
fibers  and  the  tensile  softening  in  the  matrix  process  zone 
near  the  crack  tip  can  be  easily  accounted  for  by  endo- 
chronic formulations. 

4)  Another  interesting  area  of  research  is  in  the  probabi- 
listic characterization  of  the  endochronic  parameters  by 
fitting  largely  scattered  data  from  different  sources. 

Then  the  statistical  mean  response  of  the  structure  and  the 
associated  reliability  coefficients  can  be  determined. 

Since  each  set  of  constants  represents  a particular  type 

of  behavior,  the  probabilistic  endochronic  approach  would 
provide  a good  basis  for  this  kind  of  analysis. 


6.3  Conclusions 


It  is  relevant  to  conclude  this  investigation  with 
a philosophical  remark  that  no  single  constitutive  model  can 
be  expected  to  be  valid  under  all  possible  circumstances.  The 
fact  that  it  may  not  be  valid  for  some  situations  should  not 
obscure  its  superiority  and  usefulness  in  other  cases. 


APPENDIX  A 

NUMERICAL  IMPLEMENTATION  TECHNIQUES 


A. 1 Stress-Strain  Model 

Stress-strain  curves  are  generated  by  the  endochronic 
formulation  described  in  Section  3.3.  The  curve  is  gener- 
ated by  a strain  control  procedure  in  which  the  stresses 
corresponding  to  a given  strain  increment  are  found.  Eqs. 
3.1a  and  3.1b  can  be  written  in  an  expanded  form  as  follows: 
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and 


Aa 
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ij 

ij 


elastic  stress  increments, 
prescribed  strain  increments, 


a = volumetric  stress. 

Eq.  A.  1 is  re-cast  to  have  known  quantities  on  the  right 
hand  side  and  unknown  quantities  on  the  left.  For  example, 
for  a uniaxial  case,  all  stresses  except  are  zero  and 
Eq.  A. 1 is  rewritten  in  the  following  form: 
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For  other  loading  cases  like  biaxial  or  triaxial,  Eq.  A.1 
can  be  rewritten  in  a similar  fashion. 

For  any  given  strain  increment  in  the  r-th  step, 

Eq.  A. 3 is  solved  iteratively.  Based  on  the  values  of  (r-l)th 
step,  dZ,  dZ' , dX , dX  ' , and  D's  are  determined  from  Eqs.  3.2, 
3.3,  3.4,  3.5  and  3.10  respectively.  Then  the  inelastic 
stress  increments  Aa7j  are  determined  from  Eqs.  A. 2 and 
the  elastic  increments  from  Eqs.  A. 3.  The  iteration  is 
continued  until  the  stresses  in  two  succeeding  steps  con- 
verge to  within  1.0  percent. 
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A. 2 Determination  of  Endochronic  Constants 

Endochronic  constants  are  determined  by  an  opti- 
mization process  using  an  IMSL  (international  Mathematical 
and  Statistical  Library)  subroutine,  "ZXSSQ",  based  on  the 
finite  difference  Levenberg-Marquardt  algorithm.  The  stress- 
strain  program  is  arranged  as  a subroutine,  "FUNC",  which 
computes  the  deviations  of  the  response  curves  from  the 
specified  experimental  data.  This  subroutine  is  called  by 
the  optimization  subroutine,  "ZXSSQ",  which  varies  the  endo- 
chronic constants  to  arrive  at  an  optimal  curve-fit.  The 
flow  chart  in  Fig.  A.1  shows  the  sequence  of  operations 
involved. 

A. 3 Flexural  Behavior  of  Beams 

At  each  step  the  curvature  at  the  center  of  the  beam 
is  incremented  by  ' dK'  and  the  corresponding  moment  incre- 
ment 'dM'  is  determined.  This  moment  is  prorated  to  other 
sections  along  the  length  and  the  corresponding  curvatures 
are  determined.  Based  on  the  curvatures  at  various  sections 
along  the  length  of  the  beam,  central  deflections  are  deter- 
mined from  the  area -moment  theorem.  The  procedure  used  is 
the  same  as  that  discussed  in  the  next  Section,  A. 4 

A. 4 Axisymmetric  Finite  Element  Analysis 

Plane  strain  axisymmetric  elements,  shown  in  Fig.  A. 2 
are  used  in  the  analysis.  The  displacement  field  for  the 
element  is  written  as  follows: 
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Fig.  A. 1 Flow  Chart 
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Fig.  A. 2 Axisymmetric  Element 
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where 
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shape  function  matrix,  the  elements  of 
which  are  given  as. 
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The  stress  and  strain  are  constitutively  related  as 
{do  + do"}  = [ D ] (de)  (A. 6) 

The  stiffness  matrix  for  the  elements  is  written  as  follows: 
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The  integral  in  Eg.  A.  7 is  evaluated  at  the  centroid  of  the 
elements , 

7 = ( ri  + rj  + rm  )/3  (A. 8) 

z = ( zi  + 2j  + zm  )/3  (A. 9) 

which  gives 

Kij  = 2 77  Bi  D Bj  * A (A.  1 0 ) 

where  A is  the  area  of  the  element. 

The  nodal  forces  due  to  a distributed  load  R , are 
written  as  follows: 

{ f± } = 2 tt  r R (A. 11  ) 

The  analysis  steps  for  each  load  increment  are  as  follows: 

1 ) For  each  load  increment  dP,  the  nodal  displacements  are 
found  from 

[ K ] {du}  = {dP}  (A. 12) 

Then  the  element  strain  and  stress  increments  are  found 
as  follows: 

{de}  = [ L ] {du}  (A. 13) 

and 

{da}  = [ D ] {de}  (A. 14) 

2)  The  principal  strains  at  the  center  of  each  element  are 

checked  with  the  crushing  strain  of  the  material ( input ) . 

If  crushing  occurs,  then  the  stress  vector,  a.  , + Aa . 

l-l  i 

for  the  element  is  converted  to  nodal  forces  as 
P2  = (bT  ( CTi-l  + Aai  ) d(vol) 


(A. 15) 
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This  force  system  is  redistributed  back  into  the  rest  of 
the  structure  by  step  4.  The  element  stiffness  and 
the  updated  stress  vector  are  then  set  to  zero. 

3)  If  the  principal  stress  is  smaller  than  the  tensile 

strength  of  the  material,  a , and  the  element  has  not 
cracked  before,  step  4 is  executed.  If,  however,  the 
greatest  principal  stress  is  greater  than  cr  , cracking 
is  assumed  to  take  place  during  the  current  iteration. 
For  the  case  of  no  previous  cracking  the  stress  vector 
is  evaluated  as 

= [ DS  ] + de)  (A. 16) 

If  the  element  has  cracked  previously,  then 

dai  = [ D"  ] {dei>  (A. 17) 

and 

°i  = ai-l  + dCTi  (A. 18) 

The  [ D" ] matrix  is  formed  from  [D]  in  the  same  manner  as 
shown  in  Eq.  4.10.  The  stress  level  that  can  exist  in 
the  element  after  cracking  is  computed  and  the  residual 
stress  vector  is  computed.  The  residual  stress  vector 
is  converted  to  nodal  forces  and  redistributed  in  the 
structure  in  step  4 as 


?2  ~ ) B ^ ai  “ ai_i  + d^i>  d(vol) 


(A. 19) 


4)  The  equivalent  nodal  force  vector  for  the  entire  struc- 
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If  all  the  elements  of  P are  less  than  or  equal  to  A , 
the  next  load  increment  is  applied  to  the  structure  and 
steps  1-4  are  repeated.  If  any  of  the  components  of  P 
are  greater  than  A , the  load  P is  applied  to  the  struc- 
ture and  the  element  stresses  and  strains,  and  the  nodal 
displacements  are  computed.  This  is  repeated  until  all 
the  components  of  P are  insignificant. 

All  the  above  numerical  schemes  have  been  implemented 
on  VAX  VMS-1  1/780  computer  systems  at  the  University  of  Flo- 
rida and  Florida  Atlantic  University. 
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